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1.  Report  Summary 


The  research  supported  by  DARPA  Contract  No.  N00014-88-K-0392  concerns 
the  investigation  of  the  physical  properties  of  the  components  of  a  new  solid  state 
microbattery.  This  battery  consists  of  a  lithium  metal  anode,  a  lithium  borate 
glass  separator,  and  a  cathode  of  the  layer  compound  InSe  which  readily 
intercalates  lithium  ions  and  which  is  also  an  electronic  conductor.  The  battery 
operates  by  the  passage  of  lithium  ions  between  the  electrodes  through  the  borate 
glass  separator  and  the  passage  of  electrons  between  the  electrodes  via  the 
external  circuit.  A  single  battery  of  this  type  produces  a  voltage  of  2.7  V,  can  be 
discharged  with  a  current  density  of  SOjuAcm  and  can  achieve  a  specific  energy 
of  2>OOWh  kg~^ .  It  is  the  ultimate  objective  of  this  project  to  produce  systems 
composed  of  many  microbatteries  for  applications  requiring  high  voltage,  high 
current  density  and/or  high  energy  density. 

In  order  to  achieve  this  objective,  it  is  necessary  to  have  a  fundamental 
understanding  of  fast  ion  transport  in  the  lithium  borate  glasses  and  of  both  fast 
ion  and  electron  transport  in  the  layer  compound.  This  fimdamental 
understanding  has  been  advanced  by  both  experimental  and  theoretical  work,  the 
experimental  work  being  carried  out  at  the  Laboratoire  de  Physique  des  Solides  in 
Paris  by  the  group  of  Professor  M.  Balkanski  of  the  Universite  Pierre  et  Marie 
Curie,  and  the  theoretical  work  being  carried  out  at  the  University  of  California, 
Irvine,  under  this  contract.  For  the  lithium  borate  glasses  the  experimental  work 
has  provided  measurements  of  ionic  conductivity  as  well  as  infrared  and  Raman 
spectra,  while  the  theoretical  work  has  provided  a  quantitative  understanding  of 
ion  migration  and  interionic  forces.  For  the  lithium-intercalated  layer  compound 
InSe,  the  experimental  studies  have  yielded  important  data  on  optical  properties 
due  to  electronic  transitions  and  on  infrared  and  Raman  spectra.  The  theoretical 
studies  have  yielded  the  electronic  band  structure  of  lithium  intercalated  InSe 
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and  detailed  information  on  interionic  forces  in  this  material.  The  results  of  the 
investigations  just  described  provide  the  basis  for  establishing  the  desired 
fundamental  understanding  of  fast  ion  and  electron  transport  in  the  microbattery 
components. 

Summaries  of  the  individual  projects  now  follow. 

A.  Dynamical  Properties  of  Fast  Ion  Conducting  Borate  Glasses 
Experimental  and  theoretical  investigations  which  lead  to  an  improved 

understanding  of  the  physics  of  ion  diffusion  in  lithium-doped  borate  glasses  have 
been  carried  out.  The  results  of  network  dynamics  and  lattice  dynamics 
calculations  of  the  atomic  vibrations  of  tmdoped  and  doped  borate  glasses,  i.e., 

B2O2,  and  52(93  respectively,  have  been  correlated  with  Raman  scattering 

and  infrared  reflectivity  measurements.  Addition  of  the  modifier  LioP  has  been 
shown  to,  cause  a  transformation  of  three-coordinated  to  four-coordinated  boron 
atoms.  The  force  constant  for  B-O  interactions  has  been  found  to  be  larger  for 
four-coordinated  boron  than  for  three-coordinated  boron.  The  introduction  of 
dopants,  e.g.,  LljSO^  or  UK  (withX  =  F,  Cl,  Br,  I)  leads  to  a  significant  increase  in 

the  ionic  conductivity. 

B .  Effect  of  Intercalated  Lithium  on  the  Electronic  S -ruc»  j,  of  Indium 
Selenide 

The  effect  of  intercalated  lithium  on  the  electronic  structure  of  the  (3-  and 
j-  polytypes  of  InSe  has  been  investigated  using  a  tight-binding  model.  The 
energy  bands  of  the  pure  polytypes  were  calculated  and  the  results  compared  with 
previous  work.  The  tight-binding  parameters  associated  with  intercalated 
lithium  atoms  were  obtained  using  tabulated  atomic  functions.  The  modifications 
of  the  energy  bands  and  Fermi  level  produced  by  the  introduction  of  one  lithium 
atom  per  unit  cell  were  calculated  for  the  lowest  potential  energy  positions  of  the 
lithium  atom  in  the  Van  der  Waals  gap  between  layers.  The  intercalation 
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induced  changes  in  the  smallest  and  next-to-smallest  direct  band  ^^aps  were 
determined.  An  interpretation  of  a  new  photoluminescence  peak  in  lithium- 
intercalated  InSe  has  been  developed,  and  the  existence  of  a  lithimn  impurity 
band  in  this  material  has  been  established. 

C.  Lattice  Vibrational  Properties  of  Lithium-Intercalated  indium  Selenide 

A  theoretical  investigation  of  the  normal  modes  of  vibration  of  pure  and 

lithirun-intercalated  InSe  has  been  carried  out.  Lattice  dynamical  models  have 
been  developed  for  both  materials.  The  force  constants  were  determined  by  fitting 
calculated  normal  mode  frequencies  to  experimental  infrared  and  Raman 
spectra.  It  has  been  found  that  localized  vibrational  modes  associated  with 
lithium-atom  vibrations  appear  in  the  lithium-intercalated  InSe.  One  localized 
mode  occurs  at  a  very  high  frequency  and  corresponds  to  a  lithium  atom  vibrating 
perpendicular  to  the  layer  plane,  while  a  second  localized  mode  occurs  at  a  much 
lower  frequency  and  corresponds  to  a  lithium  atom  vibrating  parallel  to  the  layer 
plane.  The  results  of  this  investigation  provide  precise  information  concerning 
the  forces  of  interaction  between  an  intercalated  lithium  atom  and  the 
neighboring  indium  and  selenium  atoms. 

D.  Molecular  Dynamics  Study  of  the  Lattice  Vibration  Contribution  to  the 

Freouencv-Dependent  Dielectric  Constant  of  Lithium  Iodide 

A  promising  approach  to  realistic  calculations  of  fast  ion  dynamics  in 
lithium  borate  glasses  involves  molecular  dynamics  whereby  one  solves  the 
classical  equations  of  motion  of  the  ions  on  a  computer  for  a  given  set  of  interionic 
potentials.  This  is  a  very  complicated  problem  for  a  disordered  system  such  as  a 
glass,  so  we  have  started  this  project  by  performing  a  molecular  dynamics 
simulation  of  a  simpler  system,  namely,  the  crystal  lithium  iodide,  Lil.  A  rigid 
ion  potential  was  used  with  parameters  fit  to  thermal  expansion,  isothermal 
compressibility,  lattice  energy,  and  the  frequency  of  the  transverse  optical  mode  at 


the  zone  center.  The  cxirrentr-cmrent' correlation  function  was  calculated  at 
T  =  200K  and  400/1^,  and  from  this  function  the  real  and  imaginary  parts  of  the 
dielectric  function  as  functions  of  frequency  were  obtained.  Anharmonic 
broadening  was  observed  at  the  higher  temperature.  The  experience  gained  from 
this  calculation  will  be  of  great  value  in  carrying  out  the  more  difficult  molecular 
dynamics  simulation  for  a  lithium  borate  glass. 


II.  Detailed  Project  Reports 


A.  Dynamical  Properties  of  Fast  Ion  Conducting  Borate  Glasses 

M.  Balkanski*,  R.  F.  Wallis,  J.  Deppe 
Department  of  Physics 
University  of  California 
Irvine,  CA  92717,  U.S.A. 

M.  Massot 

Leiboratoirc  de  Physiques  des  Bolides 
Universite  Pierre  et  Marie  Curie 
4,  place  Jussieu,  75252  Paris  Cedex  05,  FRANCE 

1.  Introduction 

Fast  ion  conducting  glasses  have  attracted  considerable 
interesr  in  recent  years  in  view  of  their  use  as  separators  in 
solid  state  batteries.  A  solid  state  cell  is  formed  of  three 
components:  an  ion  source,  usually  Li;  a  dielectric  insulator  but 
good  ion  conductor  such  as  the  lithiiam  borate  glass,  for  example, 
which  acts  as  a  separator  between  the  anode,  ion  source  and  the 
cathode  or  electron  exchanger;  and  an  insertion  compound  acting  as 
an  electron  exchanger.  The  requirements  of  the  separator  are  to 
be  shape  adaptable  in  order  to  assure  good  contact  with  the  ion 
source  and  the  exchanger  and  to  have  a  good  ionic  conductivity. 
Anocher  important  characteristic  is  that  the  separator  should  be 
as  perfect  a  dielecrric  as  possible  in  order  to  prevent  electronic 
leakage  and  assure  a  long  shelf  lifetime  of  the  batteries.  The 
borate  glasses  are  one  of  the  best  materials  among  those  known 


today  to  respond  to  all  these  reqiiirs2'*snts . 

The  boron  oxide  glass  is  a  very  good  insulator  and 

constitutes  the  glass  former.  Upon  the  addition  of  a  modifier 

such  as  Li20,  for  example,  the  structure  of  the  glass  is 

significantly  changed.  In  addition,  some  fraction  of  the  Li20 

dissociates  releasing  free  lithium  ions  and  thus  considerably 

increasing  the  ion  conductivity  of  the  glass.  A  further  increase 

of  the  conductivity  is  achieved  by  adding  a  dopant  containing  the 

mobile  ion,  Li  in  this  case.  The  effect  of  the  dopant  upon  the 

concentration  of  free  ions  available  for  the  conduction  process 

has  been  investigated  recently  [1].  Spectroscopic  investigations 

on  the  structure  of  borate  glasses  with  different  modifier 
concentrations  and  different  dopants  have  also  been  extensively 

developed  and  a  summary  is  given  in  a  recent  publication  [2]. 

The  structural  modifications  conditioning  the  fast  ion 
conduction  in  borate  glasses  are  inferred  mainly  from  spectro¬ 
scopic  data  which  reflects  the  dynamical  properties  of  the  glass 
matrix.  It  was  therefore  important  to  develop  theoretical  studies 
on  the  vibrations  of  the  glass  in  parallel  with  the  experimental 
investigations.  We  shall  present  here  simultaneously  the 
theoretical  and  experimental  results  on  the  vibrational  properties 
cf  the  borate  glass  with  the  aim  of  reaching  some  conclusions  on 
the  conditions  for  increasing  the  ionic  conductivity  of  the 
glass.  In  successive  chapters  we  shall  deal  with  the  vibrations 
of  pure  borate  glass.  Chapter  2,  and  with  the  modifications 
occurring  when  different  concentrations  of  Li20  are  added  to  the 
borate  glass.  Chapter  3.  Chapter  4  will  be  devoted  to  the  effect 
of  different  dopants,  and  in  Chapter  5  we  shall  discuss  the  effect 
of  the  glass  structure  on  the  lithium  mobility  and  present  a 


statistical  mechanical  theory  of  ionic  conductivity.  In  chapter  6 
a  double  well  model  of  frequency-dependent  conductivity  will  be 
discussed. 

2 .  Vibrations  of  Boron  Oxide  Glass 

The  structure  of  the  boron  oxide  glass  has  been  studied  by 

-'~tay  {21,  Raman  [4],  infrared  [5],  nuclear  magnetic  resonance  [6] 

and  neutron  diffraction  [7]  measurements.  According  to  these 

observations,  boron  oxide  glass  is  formed  of  an  infinite  network 

of  boroxol  groups  B3O3.  The  six  atoms  inside  the  croup  form  a 

planar  ring  and  the  B-O-B  angle  is  120°.  These  boroxol  hexagons 
are  interconnected  by  an  oxygen  bridge  -O-  and  the  external  B-O-B 

angle  i ,  less  accurately  determined,  is  about  130°.  The  relative 

geometry  of  two  successive  boroxol  planes  is  not  well  known. 

One  of  the  characteristic  features  of  the  vibrational  spectra 
of  V-B2O3  observed  by  Raman  scattering  [2]  and  shown  in  rig.  1  is 
the  presence  of  a  sharp  peak  at  800  cm“^  which  is  strong  and  well 
polarized  [31-  Usually  the  vibrational  spectra  of  amorphous  solids 
are  expected  to  consist  of  broad  bands  (see  for  example  ref.  9). 
This  is  because  the  localized  molecular  mode  responsible  for  a 
sharp  peak  couples  with  all  other  atoms  and  broadens  when  a  net¬ 
work  is  formed.  For  example,  Galeener  et  al.  [8]  have  calculated 
vibrations  of  V-B2O3  using  the  Bcthe  lattice  approximation 
(described  in  detail  below)  and  obtained  only  broad  bands  due  to 
extended  modes. 

There  have  been  a  few  attempts  to  understand  theoretically 
the  possible  existence  of  sharp  features  in  the  vibrational  spec¬ 
tra.  Kristiansen  and  Krogh-Moe  [10]  calculated  the  vibrational 
frequencies  of  the  boroxol  molecule  with  atoms  of  variable  mass 


Wavenumber  (cm-^) 

Fic.  1.  Raman  data  and  infrared  data  (inset)  for  pure  B2O3.  The  sharp 
ealc  a"  308  c-~~  is  aurributed  to  the  breathing  node  of  the  six- 
leiticered  bcroxol  rino,  i/_,  in  Fia.  2. 


attached  to  the  outer  oxygen  atoms  and  found  that  the  node  at  800 
cm"*  does  not  vary  much  when  the  mass  is  changed.  Galeener  et  al. 
i'll]  explain  the  decoupling  of  rings  fron  the  rest  of  the  network 
by  a  somewhat  fortuitous  cancellation  of  central  and  noncentrai 
forces . 

Galeener  and  Thorpe  [12]  have  also  discussed  the  introduction 
of  6-fold  borcxol  rings  into  a  continuous  random  network  model  of 
V-B2O3.  They  conclude  that  the  800  cm~”  mode  consists  of  a 
breathing  motion  of  the  intra-ring  oxygen,  and  therefore  does  not 
couple  to  the  outside  network. 

Kanehisa  and  Elliott  [13]  have  demonstrated  this  decoupling 
explicitly  using  a  cluster  Hethe-lattice  approximation  as  outlined 
in  the  following  section. 

2.1.  The  Sethe  lattice  approximation 

M.  A.  Kanehisa  and  R-  J.  Elliott  [13]  have  used  an  extension 
of  the  Bethe  lattice  approximation  to  calculate  a  partial  vibra¬ 
tional  density  of  states  of  the  boroxol  network.  They  show  by 
explicit  calculation  that  the  800  cm~~  mode  of  the  boroxol 
group  remains  sharp  even  when  an  infinite  network  is  formed.  This 
calculation  confirms  the  earlier  assumption  [4,51  that  this  mode 
involves  mainly  the  "breathing"  motion  of  oxygen  atoms  inside  the 
ring  while  the  boroxol-boroxol  coupling  on  network  formation 
involves  mainly  motion  of  the  boron  atoms  and  the  external 
(binding)  oxygen  atoms. 

Using  the  Eethe  lattice  approximation  Kanehisa  and  Elliott 
were  able  to  solve  the  equations  of  notion  for  an  infinite  system 
having  definite  local  symmetry  without  long-range  order-  The 
Bethe  lattice  considered  is  pictured  in  Fig.  3.  The  mathematical 
details  and  the  methods  used  to  treat  a  boroxol  unit  embedded  in 


(a) 


rig.  3.  (a)  Generic  threefold  Bethe  lattice;  (b)  boron 
ixide  glass  modelled  as  a  Bothe  lattice  of  boroxol  units  with 
cr:dging  oxygen  atoms. 
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the  Bethe  lattice  are  published  elsewhere  [13]. 

The  vibrational  spectrum  of  an  isolated  boroxol  molecule  is 
calculated  first.  By  adjusting  the  mode  frequency  for  the 
isolated  boroxol  molecule. to  the  experimental  value  of  800  cm”^, 
the  force  constant  q  is  determined  to  be  9.4x10^  dyn  cm~^.  The 
vibrational  density  of  states  calculated  for  the  network  is  shown 

in  Fig.  4.  A  width  7  of  100  cm“^  is  added  in  order  to  account  for 
the  bond  angle  fluctuation  around  the  ideal  values.  The  calcu¬ 
lated  vibrational  density  of  states  spectrum  compares  reasonably 
well  with  the  experimental  Raman  spectrum,  shown  in  Fig.  1.  There 
is  a  gap  at  1200-1500  cm“^  separating  the  low  and  high  frequency 
bands  in  the  calculated  spectrum  as  in  the  experimental  spectrum. 
The  higher  frequency  bands  are  mainly  due  to  boron  atom  motion 
while  those  below  the  g^ip  are  due  mainly  to  oxygen  atom  motion. 

The  most  remarkable  feature  in  this  calculation  is  the  fact  that 
the  mode  at  800  cm“^,  remains  sharp  without  any  significant 
broadening  on  network  formation.  This  is  probably  the  first 
explicit  demonstration  of  the  vibrational  decoupling  of  a 
molecular  mode  from  the  rest  of  the  network. 

3 .  Effect  of  the  addition  of  modifier  Li20  to  v-BqO-^  :  The  system 

B^Q^  —  xliX^Q 

3.1.  Experimental  results 

The  Raman  spectra  of  the  binary  glass  B203-xLi20  are  shown  in 
Fig.  5.  At  low  Li20  concentration,  x  =  0.1,  the  800  cm“^  band 
corresponding  to  the  breathing  oxygen  motion  in  the  boroxol  ring, 
the  1^2  mode  trigonal  deformation  of  the  ring  shown  in  Fig.  2,  is 
still  narrow  and  strong.  A  new  band  appears  at  780  cm“^  which 
corresponds  to  the  trigonal  deformation  of  six-membered  borate 
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Fig.  4.  Vibrational  density  of  states  of  the  boroxoi  netvor 
due  to  the  internal  oxygen  atoms. 


rings  with  one  or  two  BO^  units  ['2]  as  shown  in  Fig.  6  for  the 

tetraborate,  triborate  or  diborate  configuration.  The  band  at  500 

cra“^  corresponds  to  the  breathing  mode  1/3  (Fig.  2)  where  the  boron 

and  the  bridging  oxygen  atoms  move  in  phase.  The  mode  at  1240 

cm"*^,  is  associated  with  the  symmetric  vibration  of  the  same 

atoms.  The  band  between  1400  and  1550  cm“^  is  associated  with  the 
stretching  vibrations  of  the  B-0  bonds  in  BO3  units.  The  three 

regions  associated  with  the  stretching  vibrations  of  the  B-0  bonds 

which  link  the  boron-oxygen  network  will  be  modified  upon  addition 

of  Li20,  which  leads  to  the  transformation  of  trihedrally 

coordinated  boron,  to  tetrahedrally  coordinated  boron,  B^^, 

according  to  the  reaction 


/B-O-B^  +  Li  O 
0  0 


O 

I 

B-O-B— O 
I 

0 


+  2Li‘*’ 


For  low  1130  concentration,  the  abundance  of  B^^  increases,  (with 
increasing  modifier  concentration)  ;  as  a  consequence,  the  intens¬ 
ity  of  the  line  at  780  cm”^  increases  while  that  at  806  cm~^ 
decreases  and  finally  vanishes  at  x  =  0.3.  For  x  larger  than  0.5, 
the  frequency  of  the  B^^  band  decreases  until  it  reaches  760  cm~“ 
in  the  metaborate  glass  x  =  1.0. 

At  high  Li20  concentration  (x  >  0.5),  the  formation  of  non¬ 
bridging  oxygen  atoms  is  inferred  from  the  relative  intensity 
increase  of  the  two  bands  at  960  and  1480  cm“^. 

Two  stages  of  modification  can  be  distinguished  when  the 
concentration  of  the  modifier  increases.  During  the  first  stage, 
the  tetrahedrally  coordinated  boron  units  increase  and  lead  to  an 
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rig.  6.  Some  typical  borate  groups. 
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increase  of  the  linkage  of  the  boron-oxygen  network  up  to  x  =  0.5. 
Above  this  concentration,  formation  of  BO3  units  with  nonbridging 
oxygen  atoms  "opens"  the  borate  units  and  consequently  the  degree 
of  linkage  of  the  network  is  decreased. 

3.2  Theory 

With  the  vibrational  frequencies  of  groups  of  atoms  in  a 
given  configuration  deduced  from  spectroscopy  data,  it  is  of 
interest  to  calculate  the  vibrational  spectrum  with  some  specific 
short-range  configuration.  Toward  this  end,  we  have  employed  the 
network-dynamics  theory  and  lattice  dynamics  in  order  to  effect  an 
indirect  determination  of  the  glass  structure.  Central  force 
network-dynamics  seems  an  ideal  method  to  investigate  the  relation 
between  local  environmental  changes  and  the  nodes  of  vibrations, 
since  the  short  range  order  can  be  easily  modified. 

3.3.  Methods  of  Calculation 
3.3.1.  Network  dynamics 

The  network-dynamics  method,  introduced  by  Sen  and  Thorpe 
[14],  considers  only  nearest-neighbor  central  forces  in  the 
continuous  random  network  under  consideration;  the  other  forces, 
for  example,  long-range  Coulomb  and  angle  bending  forces,  are  much 
smaller  in  many  cases.  This  approximation  might  be  expected  to 
hold  particularly  well  in  the  optical  frequency  range  with  v;hich 
we  will  be  concernea. 

Thorpe  and  Galeener  [15]  have  put  the  central-force  model  of 
Sen  and  Thorpe  in  a  mere  general  framework  using  a  Lagrangian 
formulation. 

In  the  above-mentioned  formulation,  the  bond  directions  at 
each  atom  are  used  to  define  a  local  coordinate  system.  The 
labels,  1,  2',  etc.,  denote  the  sites,  or  atoms,  in  the  network 
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and  a,  a',  etc.,  denote  the  bonds  or  directions  at  each  site.  The 
displacement  of  the  atom  along  the  a  bond  is  written  as 

q^(i).  With  this  notation,  q^(i),q^,(i)  refer  to  displacements 
at  the  same  site,  but  in  different  directions,  and  q. (i),q. (i') 
refer  to  dispiacemen<'r3  of  two  nearest  neighbors  along  the  bond 
joining  them.  In  vase,  the  potential  energy  of  the  network 

takes  a  particul.-.r* v  form; 

V  =  |  V  [q.  (i)+q.(i')]^  ,  (3.1) 

where  a  is  the  nearest  neighbor  force  constant.  The  angular 
brackets  under  the  summation  indicate  that  each  bond  is  to  be 
counted  only  once. 

In  the  case  of  fourfold-coordinated  (or  planar  threefold- 
coordinated)  atoms,  the  form  an  dvercomplete  set,  so  a 

constraint  must  be  built  into  the  the  Lagrangian  to  ensure  that 
2^q^(^)=0  for  each  atom  Z.  One  can  accomplish  this,  by  adding  to 
the  potential  a  term 

I  I  [Z  q^(^)]^  (3-3V 

and  then  taking  the  limit  A  -  =>  at  the  end  of  the  calculation  to 
ensure  that  all  solutions  are  consistent  with  the  constraint 
r^q^(2)=0. 

The  expression  for  the  kinetic  energies  depends  on  the 

particular  local  environment  of  the  atom,  i.e.,  whether  it  has 

two,  three,  or  four  nearest  neighbors.  To  find  this  expression, 

we  consider  the  transformation  between  a  fixea  x,y,z  coordinate 
system  and  the  local  bond  direction  coordinate  system  as  shown  in 

Fig.  7.  WUth  these  expressions  in  hand,  it  is  straightforward  to 

generate  the  Lagrangian  for  any  particular  network  under  consideration. 
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Possible  ccordinane  sysrem  for  a  cvofold-coorair.area 
•ixes  defir.e  tLia  lonai  bend  directions.  ?  is  the 
e.  The  y-axis  bisects  the  angle  between  the  q  axes. 


=  T({q.))  -  V({q^))  , 


(3.5) 


where  the  curly  braces,  {...),  indicate  the  set  of  q^  for  all 
atoms  i-  The  equations  of  motion  are  then  given  by 


oL. 

dt  .• 
oq. 


3L 

3q. 


=  0 


As  discussed  in  the  previous  section,  we  have  seen  that  the 
addition  of  Li20  to  B2O3  is  accompanied  by  the  appearance  of 
groups  containing  fourfold-coordinated  borons.  In  addition,  the 
concentration  of  these  structures  increases  proportionally  with 
the  concentration  of  the  modifier  Li20  up  to  a  critical  value 
Xj,.  '.ve  shall  consider  a  network  containing  both  threefold-  and 
fouriold-coordiiiated  borons.  Indeed,  for  the  composition 
32O3-O . 5Li90,  one  can  envision  a  network  of  equal  numbers  of 
threefold-  and  fourfold-coordinated  borons.  Here  we  consider  one 
such  network,  with  each  thre-sf old-coordinated  boron  surrounded  by 
fourfold-coordinated  borons  and  vice  versa,  shown  in  Fig.  8 
Using  the  results  from  the  last  section  for  twofold-, 
threefold-  and  fourfold-coordinated  atoms,  the  Lagrangian  for  the 
network  drawn  in  Fig.  8  takes  the  form 


-2, A 


^  i 


2  A 


.3  I  \  1  2 


-I  I  fq.(^)  +  -  I  I  [y^(^)  +  ' 

where  a  =  (1-cos^)  ^  ,  a  +  2b  =  (l+cos^)”^  ,  dhci 
where  we  use  q^(i)  and  force  constant  a  for  the  fourfold- 
coordinated  borons,  y^C-?)  and  force  constant  0  for  the  threefold- 
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Fig.  B.  Network  considered  as  a  possible  srrucrure  of 
B^C; -0 .  .  :<i  IS  specified  as  being  in  the  direction  of  th 

chreefoid-coordinated  doron,  (1),  and  X2  toward  the  fourfold 
cocrdir.atec  ooron,  .  c  =  bond  label  and  therefore  direc¬ 

tor. .  've  ta.ke  2  =  force  constant  of  fourfold-coordinated  5-0 
bores,  ana  si  thar  of  threefold-coordinated  3-0  bonds. 
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coordinated  borons,  oxygen  in  the  directions  of 

q(i),y(i),  respectively,  and  again  M^oron  =  “oxygen  =  “* 
equations  of  motion  and  their  solution  are  presented  in  a  separate 
publication  [16]. 

The  central-force  eigenfrequencies  of  the  network  are 
determined  by  the  relation 

-  ma'u^)  (a^  -  I  -J)  =  (a'O^  / 

where  we  have  defined  the  frequency-dependent  fores  constants 


Mu  adCOS0 


( Mu  ^ )  ^  -  ( a +;9 )  Mu  ^ -fa  <3  S  i  n  ^  5 


a^(Mu^-5sin^3) 

Ct  .  CC^  -  '  '  '  '  ^ 

^  (Mu^)^-(a+^)Mu^4-a^sin‘'5 


5  ^  5+  _ 3^(Mu^-asin^g) 

^  (Mu^)^-  (a+/9  )Mu^^a^sin^5 


This  equation  determines  the  allowed  central- force  frequency  bands 
as  a  function  of  the  force  constants  q  and  yS,  the  masses,  and  the 
bridging  angle  9  . 

3.3,2.  Lattice  Dynamics 

The  iattice-dvnamics  calculations  were  carried  out  as  an 


independent  check  on  uhe  nerwork-dynamics  results.  As  such,  chey 
were  done  in  the  nearest-neighbor  cenrral-force  approximation. 

No  disagreement  was  found  between  the  lattice-dynamics  and 
network-dynamics  calculations.  The  uwo  methods  are  quite  ccmple- 
nentary  i.n  examining  the  effect  of  varying  various  parameters. 
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such  as  force  constants,  intermol-ecular  bridging  angles,  or 
configurations  on  the  vibrational  spectrum  of  the  solid. 

3.4.  Results  and  Discussion 

When  a  modifier,  Li20  for  example,  is  added  to  the  glass 
former  B2O3 ,  the  modes  involving  the  B-0  bonds  which  link  the 
boron  oxygen  network  are  modified.  This  can  be  understood  as 
being  due  to  the  transformation  of  threefold-  to  fourfold- 
coordinated  boron. 

As  can  be  seen  in  the  experimental  data  for  B2O3-0.lLi2O  from 
Fig.  5,  the  vibrational  spectra  are  significantly  modified  by  the 
addition  of  lithium  oxide  [18].  A  new  band  appears  at  776  cm~“  in 
the  Raman  spectrum,  and  a  set  of  broad  features  moves  up  in  the 
infrared  soectrum,  all  due  to  the  oresence  of  50^  grouos.  At  the 
same  time,  the  intensities  of  the  808-  and  7’fi0-cni“^  peaks  in  the 
Raman  spectrum  and  the  136p-cn“^  peak  in  the  infrared  spectrum  are 
seen  to  decrease  in  intensity,  relative  to  the  other  peaks,  with 
increasing  Li20  concentration.  The  new  bands  observed  in  the 
infrared  spectrum  at  800  and  1030  cm“^  are  attributed  to  BO^  bond 
stretching  in  tetraborate  groups  [ 19 ] . 

In  Fig.  9,  the  results  of  network-  and  lattice-dynamics  cal¬ 
culations  [16]  are  presented  for  p/oc  =  0.8.  Note  that  the 
general  features  in  the  spectra  such  as  the  narrowing  of  the 
middle  frequency  gap  as  compared  to  the  data  of  pure  B2O3,  are 
found  in  the  network-dynamics  results.  The  widening  of  the  high- 
frequency  gap,  presumably  due  to  the  difficulty  of  the  modes 
associated  with  BO3  groups  to  hybridize  with  those  of  the  stiffer 
BO4  groups,  is  also  found. 
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If  we  compare  the  results  in ‘Figs.  9  and  10  with  the  Raman 
spectra  shown  in  Figs.  1  and  5,  the  agreement  is  quite  striking. 
The  results  for^=  0.9  feature  a  large  frequency  gap  which  is  in 
excellent  agreement  with  the  experimental  data  for  x  =  0.  The 
j3/cC=  0.8  calculation  results  in  more  extended  frequency  bands  which 
represents  very  closely  the  spectra  displayed  in  Fig.  'i  for 

appearance  of  modes  at  960  and  1120  cm“^  in  the 
formerly  featureless  mid- frequency  band  is  "allowed"  by  the  broad¬ 
ening  of  the  bands  predicted  by  the  network-dynamics  results.  The 
peaks  at  505  and  672  cm~“  in  the  Raman  spectra  are  basically  un¬ 
affected  by  the  addition  of  the  modifier,  and  as  expected  the  low- 
and  high-frequency  band  edges  of  the  first  band  are  unchanged. 

It  is  interesting  to  note  that  such  simple  calculations 
should  yield  results  representing  well  the  tre'^ds  in  the  spectral 
distribution  of  different  vibrational  modes  in  the  binary  glass 
system  as  a  function  of  modifier  concentration.  This  result 
suggests  that  the  force  constants  used  are  reasonable,  and  again 
underscores  the  importance  of  local  order  in  understanding  the 
vibrational  response  of  covalently  bonded  glass. 

4 .  The  Effect  of  Dopants 

4.1.  Structural  modifications  induced  by  the  Dopant 

Two  types  of  dopants  which  introduce  additional  free  Li  ions 
and  increase  the  ionic  conductivity  of  borate  glasses  are  exam¬ 
ined:  the  lithium  halides  LiX  (X  =  F,Cl,Br,I),  and  lithium  sulfate 
Li2SO..  A  general  observation  is  that  the  dope  .its  do  not 
significantly  modify  the  glass  structure  which  simply  localizes 
the  anion  leaving  as  the  only  mobile  species  the  small  lithium 
cation.  The  larger  the  anion,  the  stronger  is  its  interaction 
with  the  host  network. 
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Wovenumoer  (cm~') 

rig.  10.  Allowed  frequency  bands  and  histogram  for 
O.SLi^O  as  a  function  of  the  intermolecular  bridging  angle 
6/ot  =  .  0.9(a)  and  0.7(b).  The  histograms  are  for  i  =  120° 


Raman  and  infrared  spectra  displayed  in  Figs.  11  and  12  show 
more  pronounced  changes  in  the  spectra  as  one  goes  from  F  to  ci, 

Br,  and  finally  I.  The  main  effect  is  a  decrease  in  the  intensity 
of  the  peak  at  808  cm“^  measuring  the  relative  abundance  of  tri- 
hedrally  coordinated  boron  atoms  and  an  increase  of  the  peak 

at  734  cm““  characteristic  of  the  vibration  of  tetrahedrally 
coordinated  boron  B^^.  Analysis  of  the  shape  of  the  BO^  bands, 

784  cm“^  in  Raman  scattering  and  850-1100  cm~^  in  infrared 
reflectivity,  points  toward  the  redistribution  of  these  units  into 
the  network.  For  glass  doped  with  LiF  they  are  included  preferen¬ 
tially  in  tetraborate  groups  [21],  and  the  fluorine  anions 
participate  in  the  network  forming  BO^F  and  BO2F2  units  [22]. 
Chlorine,  bromine  and  iodine  anions  enter  the  network  in  inter¬ 
stitial  positions  [18].  Introduction  of  the  dopaiit.  t>alc  results  in 
the  breaking  of  B-O-B  links  and  the  formation  of  BO^  units.  Their 
distribution  in  the  network  is  conditioned  by  electrostatic  -inter¬ 
actions  between  the  negative  anion  and  the  BO^  units.  Modifica¬ 
tion  of  the  boron-  oxygen  network  conformation  results  in  a  less 
linked  network  facilitating  free  ion  diffusion. 

In  the  ternary  glass  B203-0.57Li20-yLiX,  in  which  the  modi¬ 
fier  concentration  is  such  that  the  glass  does  not  contain  boroxol 
rings,  the  modification  due  to  the  introduction  of  dopant  salts  is 
characterized  in  the  Raman  spectrum  by  the  appearance  of  a  shoul¬ 
der  at  about  720  cm”^  on  the  770  cm”^  band  due  to  the  vibration  of 
four-coordinated  boron  atoms.  For  a  given  halide  anion  this  new 
line  shifts  toward  lower  frequencies  and  its  intensity  increases 
with  increasing  dopant  salt  concentration,  y,  as  shown  in  Figure 
13.  At  the  same  time  the  frequency  of  the  band  at  520  cm“^  due  to 
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.g.  12.  Mid-infrared  reflecrivity  of  the  ternary  halogeno- 
alassas  B-,Ot -0 . 2Li->0"'’'LiX. 
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the  in-phase  motion  of  the  bridging  oxygen  and  the  boron  atom  in¬ 
creases.  The  stiffening  of  this  mode,  which  is  higher  the  larger 
the  anion  size,  corresponds  to  a  compression  of  the  B-0-3  band. 

There  is  also  an  increase  of  intensity  of  the  bands  at  960 
and  1480  cm'*  attributed  to  the  vibrations  of  groups  .with  non¬ 
bridging  oxygen  atoms.  Sharpening  of  these  bands  and  a  shift 
toward  higher  frequencies  is  observed  when  going  from  LiCl  to  Lil 
with  the  same  concentration. 

Infrared  reflectivity  data  corroborate  the  general  findings 
from  Raman  scattering  as  to  the  effect  of  the  addition  of  a  dopant 
salt  to  the  borate  glass.  A  new  reflectivity  band  appears  at  1245 
cm“^  when  a  dopant  is  added,  as  seen  in  Figure  14,  which  shifts 
toward  lower  frequencies  when  changing  the  dopant  salt  from  LiF 
(1245  cm”^)  to  Lil  (1220  cm'^) .  This  is  consistent  with  the 
formation  of  BO3  units  with  nonbridging  oxygen  atoms.  This  result 
points  toward  a  decrease  of  the  linkage  of  the  network.  When  a 
doping  salt  is  added,  the  network  is  expanded.  The  BO^  groups  are 
transformed  into  BO3  triangles  and  BO3  with  nonbridging  oxygen 
atoms  which  can  be  connected  into  chains,  a  configuration  which  is 
consistent  with  the  softening  of  the  mode  and  which  is  favored 
when  the  size  of  the  anion  is  increased. 

These  modifications  of  the  boron-oxygen  network  of  the 
ternary  glasses  have  shown  that  in  all  the  cases  the  lithium 
halide  addition  results  in  an  expansion  of  the  network,  the  effect 
being  most  pronounced  with  Lil.  These  results  agree  with  density 
and  glass  transition  temperature  measurements. 

These  observations  are  of  prime  importance  for  the  under¬ 
standing  of  t.he  ion  transport  properties  of  the  borate  glasses. 
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The  BO^  groups  may  be  considered  as  negative  ions  with  a  large 
ionic  radius  which  provide  binding  sites  with  small  binding  energy 
for  the  mobile  lithium  ions.  At  low  modifier  content  the  increase 
of  ionic  conductivity  comes  from  the  increase  of  the  number  of  BO. 
units  and  of  the  number  of  free  lithium  ions.  But  for  higher 
oxide  content  there  is  competition  between  two  competing  mechan¬ 
isms:  an  increase  in  the  number  of  free  lithium  ions  provided  by 
the  dissociation  of  the  dopant  salt  and  conversion  of  the  BO. 

■*T 

units  into  BO3  units  with  nonbridging  negatively  charged  oxygen 
atoms  which  trap  the  lithium  ions.  In  the  range  of  large  Li20 
modifier  content  the  conductivity  is  less  sensitive  to  the 
addition  of  the  doping  salt. 

4.2.  Effect  of  the  dopant  Li2S04;  B203-xLi20-yLi2SO. 

Comparison  of  the  Raman  scattering  spectra  of  the  binary 
and  ternary  B203-0.7Li20-0.42Li2S04  systems,  dis¬ 
played  in  Figure  15 ,  shows  that • the  spectrum  of  the  Li2S04  doped 
glass  consists  of  the  superposition  of  the  spectra  of  B2O3-0.7Li2O 
and  Li2S0..  The  peaks  appearing  at  456,  644,  1004  and  1100  cm"^ 
correspond  respectively  to  i/^,  v  and  1/ ^  vibrations  of  the 

sulfate  ion  [23]. 

Analysis  of  the  infrared  reflectivity  spectra  presented  in 
Figure  16  reveals  the  disappearance  of  the  band  at  1235  cm~“  due 
to  rotations  in  the  trigonal  BO3  units  with  nonbridging  oxygen 
atoms.  Sulfate  addition  transforms  these  groups  into  BO^  units 
included  in  the  tetraborate  or  diborate  groups  [24].  Foirmation  of 
diborate  groups  preferentially  indicates  that  there  is  an 
expansion  of  the  network  to  accommodate  the  sulfate  ions  with  a 
decrease  of  the  degree  of  linkage  of  the  network. 
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g.  15.  Raman  spectra  of  the  ternary  suifoborate  glasses 
7Li-,0-'/Li2*Q4  versus  the  lithium  sulfate  content. 
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rig.  16.  Mid-i 
ate  glasses  32^2" 
tent.  The  cur'/e 
the  ratit  cf  the 
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irared  ref  ]  ectiviry  of  the  ternary  sulfo- 
7Li20-y  Li2SO^  versus  the  iithiua  sulfate 
-.eluded  i.n  the  inset  shows  the  variation 
“.tensity  of  the  30^  band  to  that  or  the  sC 


d  with  sulfate  concentration. 


Ion  conduction  in  lithium  borate  alasses. 


5.1.  Introduction 


We  have  seen  in  the  previous  chapter  that  in  this  superionic 

conductor  the  dopant  salt  LiCl  or  Li2S04  is  dissolved  in  a  solid 

matrix  consisting  of  the  lithium  borate  glass  B203-xLi20.  The 

salt  dissociates  giving  Li'*'  ions  and  Cl“  or  LiS04-  20  2- 

4 

ions.  The  Li'*’  ions  are  highly  mobile  and  significantly  increase 
the  electrical  conductivity  of  the  system.  The  situation  is  quite 

analogous  to  the  electrolytic  solution  where  a  salt  is  dissolved 
in  a  liquid  solvent  such  as  water.  There  is  one  significant 

difference,  however,  and  that  is  that  water  has  a  dielectric 
constant  of  about  80  whereas  borate  glasses  have  a  dielectric 


constant  of  about  8.  This  means  that  at  room  temperature  a 
typical  salt  such  as  LiCl  will  be  nearly  fully  dissociated  in 
aqueous  solution  but  will  be  only  weakly  dissociated  in  the  borate 
glass.  Furthermore,  in  the  latter  case,  the  concentration  of 
lithium  ions  will  increase  rapidly  as  the  temperature  is  increased 
above  room  temperature.  Consequently,  the  electrical  conductivity 
of  the  doped  borate  glasses  should  increase  rapidly  with 
increasing  temperature  as  is  well  demonstrated  experimentally. 

The  electrical  properties  of  these  glasses  have  been  studied 
by  the  method  of  complex  impedance  in  the  range  from  20  Hz  to  1 
MHz  at  temperatures  ranging  between  20  and  400°C.  The  activation 
energy  was  calculated  from  the  conductivity  a  according  to  the 
formula 
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The  values  obtained  lie  in  the  range  0.46-0.72  aV.  At  low 
concentrations  x  and  y,  the  value  of  the  activation  energy 
doorcases  rapidly  as  either  concentration  is  increased  but  changes 
only  slowly  at  high  concentrations. 

5.2,  Calculation  of  mobile  Li"^  ion  concentration  r251 
The  system  that  we  are  considering  can  be  compared  with  that 
of  an  alkali  halide  such  as  NaCl  doped  with  CaCl2  CdCl2^  which 
has' bean  discussed  previously  [26].  The  Ca^"*"  or  Cd^"^  ions 
substitute  for  Na'*’  ions  in  the  lattice  and,  to  presexrve  electro¬ 
neutrality,  positive  ion  vacancies  must  be  introduced  simul¬ 
taneously  with  a  concentration  equal  to  that  of  Ca^"*"  or  Cd^'*’. 
Positive  ion  vacancies  behave  as  negative  charges  which,  whan 
dissociated  from  the  Ca^"^  or  Cd^"^  ions,  are  mobile  and  can 
>-wutribute  to  the  electrical  conductivity.  The  overall  activation 
energy  for  the  conductivity  is  the  sum  of  a  part  connected  with 
the  dissociation  of  a  positive  ion  vacancy  from  a  divalent  cation 
and  another  part  connected  with  the  diffusion  of  a  positive  ion 
vacancy  under  the  influence  of  an  applied  electric  field.  In  our 
system  the  mobile  lithium  ions  are  analogous  to  the  positive  ion 
vacancies.  The  concentration  of  lithium  ions  as  a  function  of 
temperature  and  dopant  concentration  is  calculated  considering  an 
ensemble  of  Li2S04  molecules  dispersed  in  a  borate  glass  that  is 
assumed  to  be  a  continuous  medium  [25].  If  the  molecules  are 
regarded  as  being  isolated,  and  we  assume  that  only  one  Li"*"  ion 
can  be  ionized  per  molecule  (these  assumptions  being  valid  in  the 
low  concentration  limit  at  temperatures  which  are  not  too  high) , 
the  conductivity  can  be  written  as 
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=  ean^^/2{yn^exp(-  5^)  -  xn'exp(-  "1^)} 


where  x  =  n'/n«r. 
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where  r^  is  the  distance  of  closest:  approach  of  an  Li"^  ion  to  an 
LiO“  ion,  ^0  distance  of  closest  approach  of  an  Li~  ion  to 

an  LiSO<“  ion. 

The  lower  limit  rg  can  fae  estimated  from  Pauling's  values 
[27]  for  the  S-0  distance  and  the  radii  for  O  and  Li~.  Adding 
these  values  gives  the  result  r^  =  2.10k.  From  Pauling's 
radii  [27]  one  obtains  the  value  2.0A  for  r^. 

Using  the  Einstein  relation  connecting  the  mobility  p  and  the 
diffusion  coefficient  D,  and  recognizing  that  both  quantities  may 
have  an  activation  energy  associated  with  activated  hopping. 


one  can  write 


eP 

^  k  T 


_  ^  “ah^ 

“  kgT  jCgT*  ■ 


one  sees  that  a  plot  of  ln{c'S'^'  vs.  l/T 
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for  X  =  0  or  y  =  0  should  yield  a  straight  line.  Even  when 
neither  x  nor  y  is  zero,  this  plot  should  still  be  very  close  to  a 
straight  line.  When  the  slope  of  the  line  is  multiplied  by  -  kg 
the  effective  activation  energy  for  the  conductivity  is 
obtained. 

In  this  case  the  parameter  n-p  is  taken  to  be  i.o  x  lo^^  cm“^ 
while  the  parameters  Dq,  E-^j.  and  c  are  varied  to  give  the  best 
least-squares  fit  to  the  experimental  data.  The  mobility  activa¬ 
tion  energy  Egj^  is  determined  from  the  y  =  0  data  (with  only  Li20 
present)  to  be  0.38  eV  and  is  kept  at  this  same  value  for  the 
y  0  data  (with*  Li2S0.  present)  .  The  parameters  Dq  and  £  are 
determined  for  each  value  of  y.  Calculations  were  carried  out 
using  both  procedures.  The  results  are  as  follows. 

3/2 

The  quantity  ln(aT  '  )  is  plotted  against  the  reciprocal  of 
the  absolute  temperature  in  Fig.  17  for  each  of  the  compositions 
y  =  0.0,  0.05,  0.10  and  0.15.  The  theoretical  curves  are  straight 
lines  which  fit  the  experimental  data  quite  well  except  for  some 
deviations  at  the  two  lowest  temperatures.  From  the  values  of  £ 
determined  by  the  fits  to  the  data,  the  values  of  the  dissociation 
activation  energies  E^^  and  E^^  were  calculated 

The  overall  activation  energy  E^  was  determined  from 

the  slopes  of  the  theoretical  curves  in  Fig.  17.  The  results  for 

^0'  ■'  ^ad'  ^ad  theoretical  and  experimental  values  of  E^ 

are  presented  in  Table  1.  One  can  see  that  the  four  quantities 

'^O '  •  '  ^ad  ^ad  ^  weak  dependence  upon  the 

composition  variable  y,  and  that  their  magnitudes  are  physically 

reasonable.  The  magnitude  of  Dq  agrees  roughly  with  the  relation 
2 

-281  Dq  -  ^ where  and  are  the  hopping  frequency  and 


TABLE  1.  Conductivity  parameters  for  B203-xLi20-yLi2SO^  without 
Debye-Hiickel  effects 


y 

“^0 

(cm^- 

s"^) 

€ 

^ad 

(eV) 

F  • 

^ad 

(eV) 

Ea 

(theo¬ 

retical) 

(eV) 

Ea 

(experi¬ 

mental) 

(eV) 

o 

o 

o 

0.14 

8.1 

— 

0.44 

1.02 

0.93 

D.05 

0.16 

8.2 

0.32 

0.44 

0.86 

0.90 

:.io 

0.22 

8.4 

0.31 

0.43 

0.83 

0.83 

0.15 

0.24 

8.5 

0.30 

0.42 

0.81 

0.82 

TABLE  2.  Conductivity  parameters 

Debye-Huckel  effects 

for  B2O3 

-xLi20-yLi2S04 

with 

Dq 

£ 

Ead 

F  * 

^ad 

Ea 

Ea 

(cm^ 

(eV) 

(eV) 

(theo- 

(experi- 

s"*) 

retical) 

mental) 

(eV) 

(eV) 

3,00  0.039 

8.1 

- 

0.44 

0.93 

0.93 

3,05  0.041 

8.2 

0.32 

0.44 

0.82 

0.90 

;.10  0.056 

8.4 

0.31 

0.43 

0.80 

0.83 

3,15  0.058 

8.5 

0.30 

0.42 

0.80 

0.82 

A1 


hopping  distance  respectively.  Taking  0.2  cm^s“^  as  being  a 
representative  value  of  Dq  from  Table  1,  and  taking  [28] 

we  find  that  =  11  A;  a  not  unreasonable  value. 

The  conductivity  a  is  plotted  against  y  in  Fig.  18(a)  for  the 
temperatures  469  and  489  K,  in  Fig.  18(b)  for  temperatures  539  and 
569  K  and  in  Fig.  18(c)  for  temperatures  of  619  and  659  K.  As  in 
Fig.  17  there  are  appreciable  deviations  between  the  theoretical 
and  the  experimental  results  for  the  two  lowest  temperatures,  but 
for  the  four  highest  temperatures  the  agreement  is  good.  One  can, 
of  course,  improve  the  agreement  for  the  two  lowest  temperatures 
by  making  «  temperature  dependent. 

The  analysis  which  has  been  presented  here  differs  signifi¬ 
cantly  from  the  standard  analysis (2S»  28)  positive 

ion  vacancies  in  the  electrical  conductivity  of  dojjed  alkali 
halide  crystals.  In  the  standard  treatment  a  positive  ion  vacancy 
is  regarded  as  being  completely  free  unless  it  occupies  one  of  the 
12  next-nearest-neighbor  sites  of  a  divalent  cation.  We  make  no 
such  assumption  in  our  treatment  and  take  into  account  the  coulomb 
interaction  at  large  separations  of  an  Li"^  ion  from  an  LiO”  or 
LiSO^”  ion. 

The  theory  which  we  have  developed  for  the  dissociation  of 
Li2S04  in  lithium  borate  glasses  B203-Li20  gives  a  good 
representation  of  the  dependence  of  the  conductivity  on  the  Li2SO^ 
content  and  on  the  temperature.  By  comparing  the  experimentally 
determined  conductivity  with  the  theoretical  results  one  obtains 
information  concerning  the  activation  energy  associated  with  the 
mobility,  the  diffusion  constant  and  the  dielectric  constant. 

A  closer  examination  of  the  experimental  data  would  require 


CONCENTRATION  OF  LijSO^  y 


CONCENTRATION  OF  Li2S0^  y 


:i=n  of  the  conducriviry  s  with  Ll^SO^ 
'a)  T  =  469  K  and  489  K,  (b)  T  =  539  and 
'.d  6  59  K. 
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some  discussion  of  the  pre-exponential  factor  related  to  the 
relaxation  time  and  of  the  effect  of  interactions  between  Li"^  ions 
on  the  conductivity.  A  detailed  examination  of  the  latter  effect 
would  require  a  more  elaborate  theory  taking  into  account  the 
Debye-Hiickel  concentration  dependence  of  the  ionic  atmosphere. 

Such  a  treatment  is  presented  in  the  next  section. 

5.3.  Effect  of  Coulomb  Interactions 

The  effect  of  Coulomb  interactions  between  free  lithium  ions 
on  the  electrical  conductivity  of  borate  glass  modified  by  Li20 
and  doped  with  Li2S04  has  been  examined  [29].  The  improved  theory 
shows  a  distinct  lowering  of  the  conductivity,  at  least  for 
reiatively  low  temperatures,  arising  from  the  Coulomb  interactions 
between  free  Li"*"  ions.  At  low  temperatures  a  given  Li"*"  ion  is 
slowed  by  its  interactions  with  its  ion  atmosphere. 

The  effect  of  the  Coulomb  interaction  between  free  Li"*"  ions 
can  be  taken  into  account  by  using  the  Debye-Hiickel  theory  in  a 
form  applied  by  Lidiard  [26]  to  a  similar  problem.  What  has  to  be 
done  is  to  replace  the  dissociation  constants  K(T)  and  K'  (T)  by 
Kq_h(T)  and  K^_jj(T)  given  by 


■'d-h*'"'  =  ^  irk-}] 


=  k  «p[k:  ■  ikiH' 


where  =  8ffn_.j,.\ .  The  quantities  r^^  minimum 

distances  of  Li”  from  LiSO^”  and  LiO”  respectively,  for  which  the 
Li”  ion  can  be  considered  to  be  free.  We  do  not  know  r,^  and  rj^ 
a  priori ,  but  we  anticipate  that  rj_>rQ  and  rj^  >  r^J^.  Their  precise 


specification  will  be  given  later.  Since  the  new  equilibrium  con¬ 
stants  involve  n^rp,  the  mass  action  law  equations  are  now  trans¬ 
cendental  equations  for  the  determination  of  npip  and  must  be 
solved  self-consistently. 

If  we  introduce  the  quantities  n  =  ng  +  ng  and  n'  =  ~  ^^3 

which  are  respectively  the  total  concentrations  of  Li2S04  and 
Li20,  we  can  rewrite 

Using  conservation  of  lithium  ions,  we  have 

_  _ n _  _ n_| _ 

This  is  the  equation  that  must  be  solved  self- 
consistently  to  give  the  total  free  Li'*'  ion  concentration  npp. 

Once  r.ppi  is  found,  we  can  determine  the  conductivity  a  using 
the  relation 

<7  = 

where  u  is  the  mobility.  In  the  presence  of  the  Coulomb  inter¬ 
action  between  the  free  Li'*’  ions,  the  mobility  for  the  non- 

interacti.ng  case  must  be  multiplied  by  a  factor  r  given  by 
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-  _  nf  ri-'  f ' 

^  n+h'  ^  n+n'  ' 


where 


f  =  1  - 


f  =  i  - 


3'(2^/^-^l)  (l+/crj_)  (2^'^^+Krp 


The  expression  for  the  conductivity  then  takes  the  form 


C7  =  e^nS-nj.^. 


Using  the  Einstein  relation  between  the  rabbility  and  the 


diffusion  coefficient  D,  we  can  write 


^ah>| 

kgT  kgTJ  ' 


where  is  the  mobility  activation  energy.  Also  contributing 

to  rhe  overall  conductivity  activation  energy  are  the 
dissociation  activation  energies  E^^jj  and  for  Li2S0.  and  Li20 

respectively 

The  values  of  the  parameters  Tq  and  r^ 

were  taken  to  be  those  given  in  section  5.2:  r^  =  2.7  a  and 
r^  =  2.0  A.  The  values  of  the  other  parameters  were  determined 
as  follows.  For  the  case  y  =  0  (no  Li2S04  present)  with  free- 
lithium-ion  interactions  neglected,  the  values  of  Dq,  i  and  E^^^ 
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were  determined  by  making  a  best  fit  of  the  calculated  conduce 
tivity  to  the  experimental  data.  In  the  case  of  £  ,  the  value 
reflects  additional  experimental  data  not  considered  in  section 
5.2.  The  value  of  obtained,  O'. 48  eV,  was  ,used  in  all  the 

subsequenr  calculations.  For  each  case  with  y  h  0  and  free-ion 
interactions  neglected,  we  redetermined  Dq  and  £  to  give  the  best 
fit  to  the  data. 

When  the  interactions  betv;een  free-  Li"^  ions  were  included, 
the  values  of  Dq,  r3_  and  were  obtained  by  a  best  fit  to  the 
data,  but  the  values  of  e  were  maintained  the  same  as  in  the  non¬ 
interacting  case.  Using  calculated  values  of  the  conductivity  a, 
the  values  of  the  activation  energy  E-^  were  found  from  the  slope 
of  the  plot  of  ln(c7T^'^^)  vs.  1/T.  The  best  values  of  ri  and  r| 
were  determined  to  be  r,^  =  3rQ  and  r^  =  ^r^^.  The  values  of  the 
remaining  parameters  are  presented  in  Table  2. 

The  results  for  the  conductivity  are  presented  in  Figs. 

19-21  for  the  various-  compositions  and  temperatures  considered. 
Except  for  the  highest  temperatures  and  Li2S0^  concentrations, 
there  is  a  small  but  distinct  lowering  of  the  conductivity  which 
arises  from  the  Coulomb  interaction  between  the  free  Li"^  ions. 
Although  the  interaction  leads  to  an  increase  in  the  free-ion 
concentration  due  to  shielding  of  the  negative  ion  porential,  rhis 
effecr  is  more  than  offset  by  a  reduction  in  the  mobility.  A 
given  ion  is  slowed  by  its  interaction  with  its  ion  atmosphere.  As 
the  temperature  increases,  the  lowering  of  the  conductivity 
decreases  and  is  hardly  evident  at  the  two  highest  temperatures. 
This  behavior  can  be  understood  on  the  basis  that  an  increase  in 
temperature  increases  the  random  thermal  motion  of  the  free  Li”^ 
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Fiq.  20.  Variation  in  the  conductivity  a  as  a  function  of 
Lin50.  concentration  v  for  T  =  539  K  and  T  =  569  K:  ,  with 

inclusion  of  the  interaction  between  free  Li""  ions,  -  -  withou: 
L.-.ciusicn  of  the  interaction  between  free  Li"^  ions. 
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ions  and  reduces  the  correlations  which  arise  from  the  Coulomb 
interaction  between  them.  These  results  are  consistent  with  the 
picture  that  the  lithixun  borate  glasses  behave  as  weak  electro¬ 
lytes  with  only  a  small  fraction  of  the  snd  Li2SO.  molecules 

ionized  at  the  temperatures  employed  in  the  experiments. 

5.4 .  Summary 

The  fast  ion  conductivity  in  borate  glasses  such  as 
B203-xLi20-yLi2S04  calculated  from  a  statistical  mechanical  theory 
based  on  the  assumption  of  a  partial  dissociation  of  the  dopanr 
salt,  Li2S0.,  and  the  modifier,  Li20,  compares  well  with  the 
experimental  results.  The  temperature  dependence  of  the  conduc¬ 
tivity  obtained  experimentally  is  well  reproduced  by  the  theory. 

In  the  theoretical  treatment  given  in  section  5.2,  the  Coulomb 
interaction  between  free  Li'*'  ions  was  not  taken  into  account. 

This  theory  already  gives  the  possibility  of  obtaining  information 
concerning  the  activation  energy  associated  with  the  mobility,  the 
diffusion  constant  and  the  dielectric  constant  by  comparing  the 
experimentally  determined  conductivity  with  the  theoretical 
result. 

An  improvement  in  our  treatment  is  achieved  by  the  calcula¬ 
tions  presented  in  section  5.3  which  show  that  a  distinct  lowering 
of  the  conductivity,  at  least  for  relatively  low  temperatures, 
arises  from  the  Coulomb  interaction  between  free  Li'*'  ions.  At  low 
temperatures  a  given  Li"*"  ion  is  slowed  by  its  interactions  with 
its  ion  atmosphere.  As  the  ion  atmosphere  becomes  more  diffuse  at 
higher  temperatures,  the  effective  slowing  of  the  motion  of  an 
individual  ion  by  the  Coulomb  interaction  with  its  neighbors 
decreases. 
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A  further  refinement  in  the  theoretical  treatment  of  fast  ion 


conduction  may  come  from  the  consideration  that  the  LiO~  and 
LiSO^“  ions  fixed  at  given  sites  in  the  glass  network  may  also 
have  the  possibility  of  dissociating  to  yield  Li’’  ions  which  will 
further  contribute  to  the  conductivity.  This  process  is  less 
probable  than  the  first  dissociation,  because  its  activation 
energy  is  significantly  larger  than  that  for  the  first  dissocia- 
rion.  At  sufficiently  high  temperatures,  the  influence  of  the 
second  dissociation  would  be  indicated  by  an  increase  in  the 
apparent  activation  energy. 

5.  Lithium  motion  in  borate  glasses 

6.1.  Introduction 

The  infrared  spectrum  of  lithium  borate  glasses 
3oOT-xLioO-vLi.>SO,.  can  be  regarded  as  consisting  of  two  parts: 
the  lew  frequency  domain  which  concerns  the  charge  carrier 
dynamics  and  the  higher  frequency  region  concerning  the  host 
dynamics . 

In  sections  2  and  3,  the  dynamical  properties  of  the  host 
network  were  examined.  We  shall  now  turn  to  free  carrier 
dynamics . 

Structural  and  dynamical  propeirties  are  generally  investi¬ 
gated  by  Raman  spectroscopy  [301  and  the  ionic  conductivity  is 
measured  by  means  of  complex  impedance  spectroscopy  [251. 

Infrared  spectroscopy  is  a  complementary  method  offering  the 
ocssibiiity  of  deducing  the  frequency-dependent  conductivity  at 
low  frequencies  and  the  conduction  ion  dynamics. 

The  infrared  reflectivity  of  3203-xLi20-yLi2S0-  increases 
considerably  at  frequencies  below  500  cm"*.  This  high  reflectiv- 
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ity  corresponds  to  the  motion  of  the  Li"**  ions  and  allows  us  to 
deduce  the  frequency  dependent  ionic  conductivity.  Analysis  based 
on  a  simple  double  well  potential  model  leads  to  a  picture  in 
which  the  characteristic  vibrational  fr-jquencies  of  the  Li"^  ion  is 
distributed  over  a  large  frequency  range  corresponding  to  differ¬ 
ent  sits  configurations.  The  extent  of  the  frequency  distribution 
is  a  function  of  the  free  ion  concentration.  A  drastic  increase 
of  the  damping  coefficient  is  observed  when  the  concentration 
increases. 

The  potential  barrier  regulating  the  motion  of  the  Li"^  ions 
between  different  site  configurations  also  increases  with 
increasing  Li"^  ion  concentration  which  corresponds  to  a  decrease 
of  residence  frequency  and  increase  of  residence  relaxation  time. 

6.2.  Description  of  the  Model 

Infrared  and  Raman  spectroscopy  [30]  as  well  as  nuclear 
magnetic  resonance  measurements  [31]  indicate  that  with  increasing 
modifier  concentration  x  the  amount  of  BO^  groups  increases  up  to 
a  maximum  of  4  0%  of  boron  atoms  in  tetrahedral  coordination  for 
X  =  0.4.  Above  this  concentration  the  amount  of  BO^  decreases. 

The  relative  abundance  of  the  Bo“  ions  in  the  vitreous  matrix 
suggests  rhar  the  small  Li"^  ion  spends  a  certain  time  oscillating 
against  that  radical.  The  motion  of  the  Li"^  ion  at  long  distances 
through  the  glass  medium  can  be  characterized  by  an  activation 
energy  £_  much  larger  than  the -double-well  potential  t  (Fig.  22). 
The  Li"^  ion  should  overcome  this  potential  when  moving  from  one 
position  in  the  neighborhood  of  BO^  to  another.  Thus,  the  dynam¬ 
ics  of  the  Li'*’  ion  could  be  schematically  treated  within  the 
framewor.k  of  a  double  potential  well  model. 
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The  equation  for  the  oscillatory  motion  of  the  Li"^  ion 


considered  as  a  Brownian  particle  can  be  written  as 


m  ji  -  m^x  +  U*(x)  =  f(t)  , 


with 


U*  (X) 


3U(X) 

3X 


and 

0,X)  =  -  2{^y]  . 

o  o 

U(x)  describes  the  double  potential  well  with  minima  situated  at 
X  =  iXQ.  The  damping  coefficient  7  and  the  random  force  f(t) 
describe  the  interaction  of  the  mobile  ions  with  the  thermal 
vibrations  of  the  network. 

Following  the  treatment  given  by  Schneider  et  al.  [32], 
frequency-dependent  conductivity  can  be  written  in  the  form 


a-(w)  =  f  r  (w  ) 


Here  f  is  the  oscillator  strength  defined  as 


f  = 


“Li 


and  r  (w)  is  the  frequency  dependent  scattering  time  defined  by 
the  expression 


-1 


/ 


r (w)  =  [Z+7+r (2) ] 


with  2  =  i-j,  and  the  memory  function 


x(z) 


-2ra  l-*oi 

’■"oLz  ^  i^IbJ  • 


Here  =  <U">/m  , 


a  =  (3<X^><U''>]  , 

b  =  y^  (l-a)/ (aw^r-7  );, 


<...>  is  the  thermal  average,  and  ^  ^  =  k^T,  The:  time  r  is 

^  6SC 

related  to  the  escape  time  by  r  =  — • 

The  adjustable  parameters  for  comparison  with  experimental 
data  are  the  oscillator  strength  f,  the  lithium  vibrational 
frequency  in  the  site  the  damping  coefficient  7,  and  the 
dimensionless  parameter  S  related  to  the  height  of  the  potential 
barrier  S  =  4A^.  With  the  frequency  dependent  conductivity  <7(w) 
in  hand  we  can  calculate  the  dielectric  function 


£  (w) 


€ 

00 


+ 


Ania  (u ) 

tj 


/ 


and  confront  the  reflectivity  which  is  obtained  experimentally 

,  _  '2 

76-1 


R  = 


/X  +  l 


6.3.  Comparison  with  Experimental  Data- 

Infrared  reflectivity  measured  with  the  Bruker  type  IFS113 
vacuum  infrared  Fourier  spectrometer  give  broad  spectra  over  the 
range  of  10  -  3000  cm”^.  In  Fig.  23  is  shown  the  reflectivity 
spectrum  of  a  glass  with  composition  B2O3-0.7Li20  taken  at  room 
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temperature.  The  dotted  line  represents  the  experimental  points 
and  the  continuous  line  is  the  calculated  curve. 

Figure  24  shows  the  frequency  dependent  conductivity  a  (w)  for 
the  same  sample,  the  real  and  imaginary  parts  of  the  dielectric 
constants  and  the  energy  loss.  The  points  represent  the  Kramers- 
Kronig  inversion  of  the  experimental  reflectivity  and  the 
continuous  line  gives  the  calculated  curves.  The  maximum  of 
Im[«(v)]  in  Fig.  24  situated  at  =  40  cm~^  determines  the 
relaxation  frequency  of  the  Li"^  ion  jumping  from  one  potential 
well  to  another  through  the  barrier  A. 

When  a  dopant,  such  as  Li2S0  4,  is  added  to  the  borate  glass 

the  ionic  conductivity  increases.  A  series  of  doped  borate  glasses 

3203-0.7Li20-yLi2S04  in  a  broad  concentration  range  of  the  dopant 

salt,  0  r  y  <  0.5  has  been  studied  with  the  same  experimental 

technique.  The  remarkable  feature  in  this  investigation  is  that 

whereas  the  frequency  of  the  maximum  of  the  broad  band  w  does  not 

o 

change  with  concentration,  one  observes  a  significant  change  of 

the  damping  coefficient  7.  The  damping  coefficient  varies  from 
—  1 

7  =  240  cm  '  when  y  =  0  to  7  =  62  cm  for  y  =  0.5.  From  the 
conductivity  cr{tj),  the  real  and  imaginary  parts  of  the  dielectric 
function,  Re(«)  and  Im(£),  and  Im(“l/£)  deduced  from  the  reflec¬ 
tivity  spectra  for  different  concentrations  y  of  the  dopant 
Li2SO^,  one  can  deduce  that  the  relaxation  frequency  is  now  20 
cm“^,  a  drasric  decrease  with  increased  doping. 

From  Fig.  25,  which  is  a  summary  of  the  experimental  results, 
one  observes  that  while  the  characteristic  frequency  remains 
constant,  the  damping  decreases  and  the  oscillator  strength  f 
increases  with  increasing  dopant  concentration  y. 
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From  Che  results  shown  in  Fig.  24  one  can  see  that  the  char¬ 
acteristic  frequency  for  the  Li'*'  ion  vibration  in  a  nondoped 
borate  glass  8202-0. 7Li20  is  widely  spread  over  a  large  frequency 
range  from  iDO-600  cm“^  but  has  a  maximum  at  =  3  89  cm~^,  and 
has  the  damping  constant  7  =  243  cm”^.  The  resonance  centered  at 
3  89  cm"-^  is  very  broad  which  is  an  indication  that  the  configura- 
cion  of  the  BO^  site  attracting  the  Li'*'  ion  is  not  unique.  This 
result  suggescs  rather  that  Li"*"  is  found  in  different  sites  in  the 
borate  glass  at  different  space  positions  with  regard  to  the  Bo7 
ion.  In  each  specific  site,  Li'*'  has  a  characteristic  vibration 
but  because  of  the  large  variety  of  site  configurations,  the 
characteristic  frequencies  are  spread  over  a  broad  frequency  range 
with  a  maximum  in  the  distribution  around  390  cm"”^.  At  low  Li'*' 
concentration  the  whole  spectrum  of  different  <=ites  is  available 
for  Li"^  occupation  and  the  distribution  of  occupied  sites  is  quite 
broad.  When  the  Li'*'  concentration  is  increased  with  the  addition 
of  the  dopant  Li2S04  into  the  borate  glass  B203-0.7Li20,  the  most 
frequent  site  configurations  become  more  readily  occupied,  the 
occupancy  density  of  closely  resembling  sites  becomes  greater,  and 
the  dispersion  in  the  distribution  decreases.  The  increasing 
Li  -Li‘  Coulomb  interaction  forces  the  Li  ion  into  the  most 
favorable  configuration.  Thus,  the  resonance  sharpens  signifi¬ 
cantly  which  is  indicated  by  a  considerable  decrease  of  the 
damping  constant,  i.e.,  an  increase  of  the  relaxation  time,  the 
time  the  Li"^  spends  in  a  site  of  given  configuration.  The  oscil¬ 
lator  strength  f  increases,  of  course,  with  the  increase  of  the 
dopant  concentration  because  the  number  of  Li"*"  ions  increases. 

The  relaxation  frequency  of  Li"*",  given  by  the  maximum 


61 


of  Iin(«)  is  =  40  cm**^  for  the  nondoped  borate  glass.  This 
gives  a  potential  barrier  a  =  21.7  meV.  For  the  doped  glasses  we 
obtain  =  20  cm~^  and  a  =  40  meV.  The  most  probable  or  the  most 
abundant  potential  wells  are  twice  as  deep  for  high  Li"^  concentra¬ 
tion.  This  is  compatible  with  a  larger  residence  time  and  smaller 
damping  constants,  i.e.,  small  escape  frequency.  From  the 
definition  of  the  oscillator  strength  one  should 
expect  a  linear  relation  between  and  f  with  a  given  slope 

(e^/mLj) .  Comparison  of  the  values  of  f  obtained  from  the  fit 
with  the  experimental  data,  f^xp'  that  calculated  with  eq. 

(6.5),  fcal'  ®  considerable  difference.  Considering  the 

values  listed  in  Table  3,  one  observes  that  f^^p  practically 
double  fcal*  This  suggests  that  e  and  mj^^  are  not  the  charge  and 
mass  of  the  free  Li"^  ion  but  are  instead  the  effective  charge  (e*) 
and  effective  mass  (m*)  as  modified  by  the  polarization  induced  in 
the  medium.  The  sites  in  which  Li’*’  vibrates  are  polarizable,  which 
must  be  taken  into  account  when  calculating  f 

The  frequency  dependent  conductivity  sum  rule  is  different  for 
a  rigid  and  for  a  polarizable  medium.  For  a  rigid  medium  we  have 
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whereas  for  a  polarizable  medium  the  sum  rule  is 
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Here  e*  is  the  Szigeti  effective  charge  and  m*  the  reduced  mass 
which  for  a  tetrahedral  oxygen  environment  of  the  Li’*'  will  be 
given,  for  example,  by 
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TABLE  3.  Values  of  lithium  concentration,  dielectric  coinstant, 
and  oscillator  strength  for  various  concentrations  y  of  the  dopant 


• 

y 

Nj^j_(cm“^) 

€ 

oo 

^exp/^cal 

0.0 

2.2  10^^ 

2.06 

.21774 

12204 

1.784 

0.1 

2.28  10^2 

2.06 

24704 

12648 

1.953 

• 

0.2 

2.32  10^2 

2.18 

27369 

12869 

2.127 

0.4 

2.38  10^2 

2.13 

29512 

13202 

2.235 

nij.. 


The  ratio  e  /e  is  given  by 


For  different  dopant  concentrations  y  we  obtain  the  values  given 
in  Table  4* 


This  result  shows  that,  in  fact,  the  polarizability  of  the 
^  3203-0.7Li20-Li2S04  network  does  not  affect  considerably  the  Li’*' 

oscillatory  motion  because  the  effective  charge  e*  is  not  very 
different  from  the  free  Li"^  charge  e. 
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TABLE  4. 


Effective  charge  of  Li’’*  for  various  concentrations  y  of 


the  dopant. 


y 

e*/e 

0 

0.81 

0.1 

0.85 

0.2 

0.86 

0.4 

0.89 

7 .  Conclusions 

After  long  discussions  in  the  literature  it  now  seems 
established  that  the  network  is  mainly  formed  by  boroxol 

rings  connected  with  B-O-B  bonds  and  BO3  triangles.  Direct 
calculation  of  the  vibration  of  boron  oxide  glass,  using  the  Bethe 
lattice  approximation,  show  that  the  boroxol  ring  vibration  at  800 
cm““  remains  sharp  even  when  an  infinite  network  is  formed. 

The  effect  of  the  addition  of  the  modifier  Li20,  was  investi¬ 
gated  experimentally  by  Raman  and  infrared  spectroscopy  and  theo¬ 
retically  by  network  dynamics  and  lattice  dynamics  calculation. 

The  modifier  was  found  to  increase  the  number  of  free  lithium  ions 
which  act  as  charge  carriers,  and  also  to  significantly  modifify 
the  host  structure  through  the  transformation  of  trihedrally 
coordinated  boron  atoms  into  tetrahedrally  coordinated  boron.  The 
concentration  dependence  shows  that  this  second  effect  goes 
through  a  maximum  for  x  =  0.4  and  then  decreases  as  a  result  of 


the  formation  of  BO3  groups  with  non-bridging  oxygen  atoms  which 
carry  a  negative  charge.  Compared  with  the  experimental  data,  the 
results  of  network  and  lattice  dynamics  calculations  seems  to  give 
an  adequate  picture  of  the  structure  and  vibrational  response  of 
lithium  doped  borate  glasses. 

Introduction  of  dopants  in  the  form  of  alkali  halide  salts 
produces  by  dissociation  an  additional  free  lithi\2m  concentration 
and  also  a  modification  of  the  host  configuration  dependent  on  the 
size  of  the  anion.  Flourine  enters  in  the  network  in  substitu¬ 
tional  positions  forming  BO3F  and  BO2F2  units.  Chlorine,  bromine 
and  iodine  enter  interstitial  positions  which  produces  breaking  of 
3-0-B  links  and  formation  of  BO^  units.  This  modification  of  the 
boron  oxygen  network  results  in  a  less  linked  network  facilitating 
the  ^ree  ion  diffusion.  The  effect  is  more  pronounced  for  Lil  a»iu 
decreases  with  the  anion  size. 

Doping  with  Li2S0.  leads  to  analogous  effects.  The  effect  of 
Li2S0^  on  the  Raman  spectra  can  be  understood  as  simply  a  super¬ 
position  of  the  vibrations  of  the  host  network  and  that  of  the  SO. 
groups.  The  modifications  of  the  host  network  with  preferential 
formation  of  diborate  groups  tends  to  produce  an  expansion  of  the 
network  in  order  to  accommodate  the  sulfate  ions  and  facilitate 
the  lithium  mobility. 

The  ion  conductivity  of  lithium  borate  glasses  has  been 
addressed  by  developing  a  statistical  mechanical  theory  of 
conductivity  based  on  the  assumption  that  the  dopant  behaves  as  a 
weak  electrolyte  which  dissociates  to  provide  free  Li"^  ions  and 
high  ionic  conductivity.  The  agreement  between  the  theoretical 
results  and  the  experimental  data  is  good  for  physically  reason- 


able  values  of  the  parameters  that  characterize  the  theory. 

Taking  into  account  the  Coulomb  interactions  between  free  Li"^ 
ions  shows  a  distinct  lowering  of  the  conductivity,  at  least  for 
relatively  low  temperatures.  This  result  arises  from  the  Coulomb 
interaction  between  free  Li"*"  ions  and  the  fact  that  at  low 
temperatures  a  given  Li"^  ion  is  slowed  by  its  interaction  with  its 
ion  atmosphere. 

Lithium  motion  in  borate  glasses  has  also  been  investigated 
by  far  infrared  spectroscopy  and  analyzed  on  the  basis  of  a  double 
well  model.  This  analysis  shows  that  the  characteristic 
vibrational  frequencies  of  the  Li"*"  ion  is  distributed  over  a  large 
frequency  range  corresponding  to  different  site  configurations. 

The  extent  of  the  frequency  distribution  is  a  function  of  the  free 
ion  t^^jicentration.  A  drastic  decrease  of  the  damping  ccefficieni. 
is  observed  when  the  concentration  increases. 

The  potential  barrier  regulating  the  motion  of  Li^  ions 
between  different  site  configurations  also  increases  with 
increasing  Li^  ion  concentration  which  corresponds  to  a  decrease 
of  residence  frequency  and  increase  of  residence  relaxation  time. 

In  summary  the  experimental  and  theoretical  results  presented 
here  gives  a  fairly  broad  and  clear  picture  of  the  main  features 
of  both  the  matrix  dynamics  and  the  lithium  free  ion  dynamics  of 
lithium  doped  borate  glasses. 
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B.  Effect  of  Intercalated  Lithium  on  the  Direct  Band 


Gaps  of  Indium  Selenide 


P.  Gomes  da  Costa,  M.  Balkanski*  and  R.  F.  Wallis 
Physics  Department,  University  of  California,  Irvine,  CA  92717 

Abstract 


The  effect  of  intercalated  lithium  on  the  and  7 -polytypes 
of  InSe  has  been  investigated  using  a  tight -binding  method.  The 
energy  bands  of  the  pure  polytypes  were  calculated  and  the 
results  compared  with  previous  work.  The  tight-binding 
parameters  associated  with  intercalated  lithium  atoms  were 
obtained  using  tabulated  atomic  functions.  The  modifications  of 
the  energy  bands  and  Fermi  level  produced  by  the  introduction  of 
one  lithium  atom  per  unit  cell  were  calculated  for  the  lowest 
potential  energy  position  of  the  lithium  atom  in  the  Van  der 
Waals  gap  between  layers.  The  intercalation  induced  changes  in 
the  smallest  and  next-to-smallest  direct  band  gaps  were  deter¬ 
mined  and  compared  with  experimental  data. 
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1 .  Introduction 


The  insertion  of  lithium  into  layered  materials  has 
attracted  significant  interest  as  a  consequence  of  the  possible 
use  of  such  materials  as  electron  exchange  electrodes  in  solid 
state  batteries.  Such  an  application  is  related  to  the  inter¬ 
calation  capacity  and  the  modifications  of  the  electronic 
properties  of  the  bulk  material.  Intercalated  InSe^^^  gives  a 
voltage  of  2.6V  against  a  lithium  anode  and  exhibits  a  change  in 
conductivity  of  up  to  two  orders  of  magnitude. 

InSe  is  a  1  ^.yered  material  consisting  of  two-dimensional 
Se-In-In-Se  layers  between  which  the  binding  is  weak  due  to  the 
Van  der  Waals  forces.  The  space  between  successive  layers  is 
known  as  a  Van  der  Waals  gap.  The  stacking  of  the  layers  is 
always  compact,  but  there  are  four  possible  stacking  arrangements 
leading  to  the  four  polytypes  /9 ,  t  ,  7 ,  and  5  .  Only  three  of 
these  polytypes  (/9,  «,  7)  have  been  observed  for  InSe.  Bridgeman 
grown  crystals  are  generally  of  the  type  7.  The  difference  in 
total  energy  between  the  various  polytypes  is  so  small  that  one 
practically  always  finds  a  high  density  of  stacking  faults  in  a 
given  sample. 

Although  most  of  the  experimental  results  that  have  been 
obtained  are  for  7 -InSe,  no  band  structure  calculations  are 
available  in  the  literature  for  this  polytype.  All  existing 
calculations  are  for  the  and  e  polytypes  which  have  a  simpler 
first  Brillouin  zone.  The  first  band  structure  calculation  for 
this  family  of  materials  was  made  for  GaSe. A  calculation  for 
InSe  within  a  tight-binding  model  has  been  carried  out  in  a  two- 
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dimensional  approximation. Extensions  to  the  three- 
dimensional  case  have  been  developed  using  the  pseudopotential 
method and  the  tight-binding  method. Including  the  spin- 
orbit  interaction  in  the  tight-binding  calculation has  a  non- 
negligible  effect  on  the  InSe  band  structure.  Pseudopotential 
calculations  for  e-InSe  are  also  available, but  without  spin- 
orbit  interaction. 

Lithium  insertion in  InSe  has  been  shown  to  affect 
both  the  electronic  and  optical  properties.  Possible  inter¬ 
stitial  sites  for  Li  atoms  in  7 -InSe  have  been  discussed  and 
their  relative  energies  determined. Lithium  diffusion 
paths in  7 -InSe  have  also  been  determined  by  ^  initio 
calculations. 

Lithium  insertion^^' in  InSe  has  been  shown  to  affect 
both  the  electronic  and  optical  properties.  Possible  inter¬ 
stitial  sites  for  Li  atoms  in  7 -InSe  have  been  discussed  and 
their  relative  energies  determined. Lithium  diffusion 
paths in  7 -InSe  have  also  been  determined  by  ^  initio 
calculations. 

The  optical  properties  near  the  band  edges  of  pure  InSe  have 
been  investigated  and  compared  with  the  calculated  band 
structure. Sharp  excitonic  peaks  are  observed  in  the 
absorption  spectrum  at  low  temperature  corresponding  to  the  three 
absorption  thresholds.  The  first  threshold  at  1.3  eV  is  related 
to  the  direct  absorption  between  the  s,  valence  band  states 
and  s  conduction  band  states.  The  second  at  2.5  eV  is  related  to 
transitions  from  Pj^,  Py  valence  band  states  to  s  conduction  band 


states,  and  the  third  at  2.9  eV  is  related  to  transitions  from 
the  spin-orbit  split-off  valence  band  to  the  s  conduction  .iid. 
More  recent  results on  pure  InSe  at  lOK  with  the  electric 
vector  of  the  radiation  perpendicular  to  the  c-axis  allow  one  to 
distinguish  the  n=l  and  n=2  excitonic  transitions  as  well  as  the 
LO  phonon  replicas  of  the  n=l  excitonic  state. 

The  effect  of  Li  insertion  on  the  interband  optical 
absorption  peaks  is  rather  weak,  but  nevertheless,  clearly 
observable.  In  Fig.  1  are  shown  the  shifts  in  frequency  of  the 
smallest  direct  band  gap  and  the  next-to-smallest  band  gap  E2 
as  functions  of  the  lithium  content.  We  see  that  the  E^^  gap 

increases  in  frequency  as  the  Li  concentration  increases,  whereas 
the  E2  gap  decreases  in  frequency. 

The  excitonic  transitions  persist  after  Li  insertion,  which 
suggests  that  all  of  the  Li  2s  electrons  do  not  transfer  to  the 
conduction  band  and  thus  transform  the  semiconducting  InSe  into  a 
metal. If  we  had  a  metallic  transition,  the  Coulomb 
interaction  between  the  electron  and  hole  of  the  exciton  would 
be  screened  and  the  excitonic  state  would  be  washed  out.  The 
persistence  of  the  excitonic  transitions  in  highly  intercalated 
InSe  suggests  that  the  Li  2s  electrons  form  a  low  mobility 
impurity  band  or  are  efficiently  trapped  into  localized  states. 

The  photoluminescence  spectrum  of  InSe  is  significantly 
modified  by  the  intercalation  of  Li-  A  new  photoluminescence 
peak  appears  at  a  photon  energy  somewhat  less  than  that  of  the 
fundamental  exciton  peak  of  pure  InSe  as  shown  in  Fig.  2.  This 
new  peak  may  be  associated  with  the  Li-2s  band  lying  in  the 


ENERGY  (eV) 

FIG.  2.  Photoluminescence  spectra  for  7-InSe  before  (a)  and 
after  (b)  lithium  intercalation  at  5  K  under  1.916  eV  excitation 
by  a  Kr"^  laser. 
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fundamental  gap  of  InSe. 

Theoretical  Development 

We  have  developed  a  tight-binding  scheme  for  calculating  the 
electronic  band  structure  of  both  /3-  and  7-InSe.  Our  procedure 
is  based  on  the  overlap-reduced  semi-empirical  tight-binding 
method  of  Doni  et  al.^^®^  who  considered  the  case  of  p-lnSe. 

The  Slater-Koster  procedure was  used  to  express  the  overlap 
and  interaction  integrals  of  the  tight  binding  method  in  terms  of 
the  basic  overlap  integrals  S j^j  and  interaction  integrals  j . 
Only  two-center  integrals  were  considered. 

For  )9-InSe  we  took  the  values  of  the  parameters  and  Vj^j 
from  the  tabulation  in  Doni  et  al.  for  the  nearest-neighbor 
In-In,  In-Se,  and  Se-Se  interactions  both  parallel  and  perpen¬ 
dicular  to  the  c-axis.  The  atomic  states  considered  are  Se  4s 
and  4p  and  In  5s  and  5p  whose  energies  were  taken  from  Doni  et 
al.  We  also  used  the  crystal  field  parameters  K  tabulated  by 
Doni  et  al.  Interactions  beyond  nearest  neighbors  were  included 
by  scaling  the  Doni  parameters  with  a  factor  d“^  where  d  is  the 
interatomic  distance  of  the  interacting  atoms.  The  interactions 
were  cut  off  at  d=8.77A. 

The  structure  of  7-InSe  is  simpler  than  that  of  0-lTiSe, 
because  the  former  has  only  four  atoms  per  unit  cell,  whereas  the 
latter  has  eight  atoms.  The  difference  in  structure  arises  from 
the  difference  in  stacking  of  successive  layers;  the  interatomic 
distances  within  a  layer  are  essentially  the  same  in  both  poly¬ 
types.  In  view  of  these  similarities,  we  have  taken  the  basic 
overlap  and  interaction  parameters  for  7-InSe  to  be  the  same  as 


those  given  by  Doni  et  al.  with  the  exception  of  the  interlayer 
parameters  Sgg_gg(ppcr)  and  Vg^^g^Cppa)  .  These  parameters 
together  with  the  crystal  field  parameters  K(In  5s) ,  K(Se  Pj^)  , 
K(Se  Py) /  and  K(Se  P2)  were  varied  in  order  to  reproduce  the 
experimental  values  of  the  smallest  band  gap  and  the  next-to- 
smallest  band  gap  £2*  The  values  thus  obtained  are  tabulated  in 
Table  I. 

Having  established  the  tight -binding  scheme  for  the  pure  fi- 
and  7-polytypes  of  InSe,  we  then  proceeded  to  generalize  the 
scheme  to  include  the  presence  of  intercalated  Li  atoms.  We 
assume  that  the  Li  atoms  occupy  the  sites  of  lowest  potential 
energy^^^)  (the  A3_2  sites)  in  the  Van  der  Waals  gap  between 
layers.  Since  intercalation  of  lithium  leads  to  only  a  very 
small  increase  in  the  interlayer  spacing,  we  have  neglected  any 
change  in  the  tight -binding  parameters  of  InSe  itself.  The 
overlap  and  interaction  parameter  for  the  interaction  of  a  Li 
atom  with  the  nearest  Se  and  In  atoms  were  calculated  using  the 
tabulated  atomic  functions  of  dementi  and  Roetti.^^®)  The 
electronic  band  structure  of  the  intercalated  InSe  was  then 
obtained  by  diagonalizing  the  Hamiltonian  matrix.  The  changes  in 
the  E3  and  E2  gaps  thus  obtained  turned  out  to  be  large.  The  Li 
parameters  were  then  varied  until  agreement  was  achieved  between 
the  theoretical  and  experimental  gaps.  The  resultant  values  for 
the  Li  parameters  are  listed  in  Table  II. 

Numerical  Results 

The  electronic  band  structure  of  pure  jT-InSe  which  we  have 
obtained  is  shown  in  Fig.  3  for  various  directions  in  the  first 
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TABLE  I.  Values  of  the  tight-binding  parameters  for  7-InSe.  The 
notation  is  that  of  Ref.  16.  V  and  J  are  in  Rydbergs. 


m 

Se-Se  interlayer  parameters 

crystal  field 

parameters 

W' 

S(ppff)  -0.010 

K(Se  4Pjj) 

0.0203 

V  (pp<r)  0.002 

K(Se  4P2) 

0.1143  1 

K(In  5s) 

-0.0373 

TABLE  II.  Values  of  the  tight-binding  parameters  for  lithium- 
intercalated  7-InSe.  V  and  K  are  in  Rydbergs. 


# 

S 

V 

K 

Li-Se 

(sstr) 

0.02636 

-0.000318 

Li-Se 

(spa) 

-0.0077 

0.0027 

• 

Li-In 

(SSa)  . 

0.004143 

-0.000313 

Li-In 

(spa) 

-0.04665 

0.000669 

• 

Li  (2s) 

—0  *  066 
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HG.  3.  Energy  bands  for  pure  y-InSe. 
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Brillouin  zone.  The  smallest  direct  gap  (Ei)  occurs  at  r  between 
the  s,  P2  valence  band  and  the  s  conduction  band.  The  next  lower 
direct  gap  (£3)  occurs  at  r  between  the  Pj^,  Py  valence  band  and 
the  s  conduction  band. 

For  the  case  of  pure  ^-InSe,  we  plot  our  results  for  the 
band  structure  in  Fig.  4,:  These  results  are  very  similar  to  those 
of  Doni  et  al.  In  particular,  both  sets  of  results  give  an 
indirect  gap  between  r  and  M  that  is  somewhat  smaller  than  the 
direct  gap  at  r . 

In  Fig.  5  we  present  the  band  structure  for  lithitua- 
intercalated  7-InSe  with  the  composition  Lig  gInSe  and  the  Li 
atoms  occupying  the  ^2-2  "the  Van  der  Waals  gap  between 

layers.  The  band  associated  with  the  Li  2s  state  is  seen  to  lie 
just  below  the  conduction  band  edge.  The  effect  of  the 
intercalated  lithium  on  the  InSe  bands  is  too  small  to  be  observ¬ 
able  on  the  scale  of  the  figure.  However,  the  effects  on  the  Et 
and  E2  gaps  are  revealed  by  the  numerical  results  given  in 
Table  Iff.,  The  excellent  agreement  between  the  experimental  and 
theoretical  values  is  a  direct  consequence  of  the  choice  of  the 
Li  parameters. 

The  occurrence  of  the  Li-2s  band  in  the  fundamental  gap  of 
InSe  provides  a  natural  explanation  for  the  new  photoluminescence 
peak  produced  by  Li  insezrtion.  Excitation  with  photons  having 
energy  greater  than  the  band  gap  creates  electron-hole  pairs.  An 
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FIG.  5.  Energy  bands  for  lithium-iniercaiated  7-IhSe  with 
the  composition  Lio.sInSe.  The  gap  between  the  Li  2s  band  and 
the  conduction  band  just  above  it  is  0.036  eV  at  F,  but  is  exag’ 
gerated  for  clarity. 


Table  III.  Values  of  the  and  E2  energy  gaps  with  and 
without  lithium  intercalation  in  7-InSe. 


pure  7-InSe 
1.3335  eV 


lithium-intercalated  7-InSe 
1.3368  eV 


2.5451 


2.539 


electron  in  the  Li-'2s  band  can  then  recombine  with  a  hole 
accompanied  by  the  emission  of  a  photon  whose  ener^  is  equal  to 
the  difference  between  the  valence  band  edge  and  the  lower  edge 
of  the  Li-2s  band. 

Conclusions 

A  tight-binding  scheme  has  been  developed  for  pure  ^-InSe 
and  7-InSe  using  the  overlap-reduced  semi-empirical  method  the 
procedure  satisfactorily  accounts  for  the  direct  energy  gaps  of 
the  two  polytypes.  The  tight -binding  scheme  has  been  generalized 
to  include  the  effect  of  intercalated  lithium  atoms.  The  shifts 
in  the  direct  gaps  and  the  appearance  of  a  new  peak  in  the  photo- 
luminescence  spectrTim  predicted  by  the  calculations  are  in 
agreement  with  experimental  optical  data. 
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ABSTRACT 


The  lattice  dynamics  of  the  7-polytype  of 
the  la-yer  compound  InSe  has  been  investigated 
using  a  model  containing-  short-range  central 
forces  and  long-range  Coulomb  interactions. 
The  normal  mode  frequencies  and  eigenvectors 
were  determined  by  diagonalizing  the 
dynamical  matrix.  The  results  have  been 
correlated  with  the  infrared  and.  Raman 
spectra  of  7-InSe.  The  effect  of 
intercalated  lithium  atoms  on  the  vibrational 
modes  has  been  investigated  both 
theoretically  and  with  the  aid  of  Raman 
spectra . 


1 .  INTRODUCTION 


In  a  number  of  solid  state  batteries  currently  under 
investigation,  layered  materials  are  used  as  electron 
exchange  electrodes,  since  they  can  accommodate  considerable 
amounts  of  intercalated  lithium.  For  such  an  application  it 
is  essential  that  the  electrode  material  have  a  sufficiently 
high  ionic  as  well  as  electronic  conductivity.  A  material 
that  satisfies  these  requirements  is  InSe  which  can 
intercalate  relatively  large  amounts  of  lithium  with  only  a 
slight  expansion  of  its  lattice  (1,2). 

InSe  is  a  layered  compound  having  complex  layers  in 
which  atomic  layers  of  Se,  In,  In,  and  Se  are  bound  together 
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by  covalent  bonds  with  some  ionic  character.  These  complex 
layers  are  themselves  bound=  to  one  another  by  weak  Van  der 
Waals  interactions  with  a  so-called  Van  der  Waals  gap 
between  layers.  There  are  four  possible  stacking 
arrangements  of  the  complex  layers  leading  to  four  polytypes 
designated  /3,e, y,  and  S  i  Since  the  experimental  results  to 
be  discussed  are  for  the  y-polytype,  we  shall  focus  our 
attention  on  this  polytype. 

A  microscopic  understanding  of  both  the  electronic  and 
ionic  conductivity  of  lithium-intercalated  y-InSe  requires  a 
knowledge  of  the  lattice  dynamics  of  this  system.  Of 
particular  importance  is  information  concerning  the  forces 
of  interaction  between  intercalated  lithium  atoms  and  nearby 
selenium  and  indium  atoms . 

In  this  paper  we  present  the  results  of  a  theoretical 
investigation  of  the  lattice  dynamics  of  y-InSe  both  with 
and  without  intercalated  lithium.  We  develop  the  lattice 
dynamical  model  in  section  2.  The  values  of  the  force 
constants  and  calculated  phonon  dispersion  curves  for  pure 
y-InSe  are  presented  in  section  3.  Similar  results  are 
presented  in  section  4  for  lithium  intercalated  y-InSe. 
Conclusions  are  presented  in  section  5. 

2.  LATTICE  DYNAMICAL  MODEL 

The  first  step  in  setting  up  the  lattice  dynamical 
model  is  to  summarize  the  pertinent  information  concerning 
the  structure  of  InSe .  The  y-polytype  belongs  to  the 
trigonal  system,  and  its  space  group  is  .  The  primitive 
translation  vectors  can  be  written  as  [3] 


{■^/3)aJ  +  c,k 


(la) 
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# 


# 
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where  is  the  lattice  constant  parallel  to  the  layers,  c^ 
is-  the  lattice  constant  perpendicular  to  the  layers,  and 
i  ,J:,k  are  unit  vectors  in  the  x,  y,  z  directions, 
respectively.  The  values  of  a^  and  c^  are  [4]  4.00  A  and 
8.44  A  respectively.- 

The  primitive  translation  vectors-  of  the  reciprocal 
lattice  can=  be  obtained  from  and  -are  found  to  be 


=  T7i[{2/Sa,)r+{\/3c,)k] 

(2a) 

=  2n[-{l  /  VSfljr  +  +  (1  /3c^k] 

(2b) 

=  27t[-{]/4ia,)T -(I7ajj  +  (l/3c,)r] 

(20 

Diagrams  of  the  stacking  of  the  layers  of  /-TnSe  and  its 
first  Brillouin  zone  are  shown  in  Fig.  1. 

Since  InSe  is  a  polar  semiconductor,  the  chemical  bonds 
which  bind  neighboring  atoms  together  have  both  covalent  and 
ionic  character.  The  interatomic  forces  should,  therefore, 
have  a  short-range  part  associated  with  the  covalent 
character  of  the  bonds  and  a  long-range  Coulomb  part 
associated  with  the  ionic  character  of  the  bonds .  We 
therefore  write  an  element  of  the  force  constant  matrix 
^afiQK;l’K')  as  the  sum  of  two  contributions, 

1^)  =  +  cD^(/«:;/'  K')  ^  ( 3 ) 
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(b) 


Fig.  1  (a)  Stacking  of  layers  and  (b)  first  Brillouin 

zone  of  y-InSe.  In  (a)  the  open  circles  are  Se  atoms  and 
the  closed  circles  are  In  atoms. 

where  is  the  contribution  from  the  long-range 

Coulomb  interactions,  is  the  contribution  from 

short-range  non-Coulomb  interactions,  the  index  pair  Ik 
designates  the  Kth  atom  in  the  Ith  unit  cell,  and  a,P  denote 
Cartesian  coordinate  components . 

For  the  non-Coulomb  interactions  we  adopted  the 
central-potential  model  which  Polian  [5)  found  satisfactory 
for  the  related  material  GaS.  With  a  central  potential  ^(r) 
dependent  only  on  che  distance  r  between  the  interacting 
atoms,  one  can  write  the  contribution  to  for 

Ik^I'k"  in  the  form 
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where 


(m> 


(5): 


x{lK;l':K’)  =  x{lK)-xirK’)  ,  (6); 

x{Ik)  is  the  position  vector  of  atom  Ik  in  the  equilibrium 
configuration,  and  primes  denote  differentiation  with 
respect  to  argument .  A  given  interaction  consequently 
involves  two  force  constant  parameters,  ^A^yir, 

which  will  henceforth  be  designated  by  A  and  B,. 
respectively . 

The  short-range  model  employed  includes  the  following 
interactions : 

1.  In-In  first  neighbors 

2.  In-Se  first  neighbors 

3 .  Se-Se  interlayer 

4.  In-In  second  neighbors 

5.  In-Se  second  neighbors 

The  corresponding  values  of  A  and  B  are  designated  by  A^^,  3,- 
with  i  =  1,2,3, 4,5. 

The  long-range  Coulomb  interactions  are  characterized 
by  the  effective  charges  of  the  In  and  Se  ions.  Since  the 
In-Se  bond  is  partly  ionic  and  partly  covalent,  we  expect 
the  magnitudes  of  the  effective  charges  to  be  less  than  the 
magnitude  of  the  electron  charge.  From  the  symmetry  and 
stoichiometry  of  /-InSe  and  electroneutrality,  we  conclude 
that  the  effective  charges  of  the  In  and  Se  ions  are  equal 
in  magnitude  and  opposite  in  sign.  The  large  anisotropy  Of 
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the  crystal  leads  us  to-  introduce  a  tensorial  effective 
charge  of  the  ;forin  [5] 

'z;  0  o' 

z*e,=  0  z^  6  e  ;  (7) 

0  0  Z2_ 

This  form,  for  the  effective  charge  accounts  for  the  fact 
that  the  polarization  produced  by  displacements  parallel  to 
the  layers  is  different  from  that  produced  .by  displacements 
.perpendicular  to  the  layers .  ' 

The  contribution  of  the  long-range  Coulomb  interactions 
to  the  force  constant  matrix  was  evaluated;  using  the  Ewald 
•technique  by  expressing  the  elements  of  this  matrix  in  terms 
of  a  Sum  over  the  direct  lattice  and  a  sum  over  the 
reciprocai  lattice v  The  elements  of  the  Coulomb 
contribution  to  the  dynamical  matrix  were  then 

■obtained  using  the  relatioii 
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where  is  the  effective  charge  compcinent  of  ion  K  given 
by  /  T]  was  taken  to  have  the  value  4fjt  lr^,  r^  is  the 
equilibrium  nearest-neighbor  separation, 


^(/;  kjS)  =  r^{K)-x{lK’% 


(10) 


ana 


JaC, 


2 


-d. 


ap 


bC 


y£ffc{C)  '^  ?  ^ 


(11) 


The  reduced  effective  charge  components  Z/^  and  satisfy 
z’s,a  =  -z]^  =  z*  with  z*  specified  by  Eq.  (7)  —  i.e.. 


•  •  •  •  • 

z,  =Zy  =  z,  and  z,  =Zj  . 


The  contribution  of  the  sum  over  the  reciprocal  lattice 
:o  is  given  by 


-—I 

'  V.  [VwA 


-k^tw 


p+kf  ^ 


where  the  reciprocal  lattice  vectors  G  are  given  by 


G-=  /Wj^i  +  ^^2  "^^^3  (4?  ) 

hnd  m^,  m^,  are  integers . 

The  full  Coulomb  contribution  to  the  dynamical  matrix 
is  given  by 

=D^^[fcvKK)  =  +  p^^(k;KK)^  (1=4 ) 

Adding  this  result  to  the  short-range  contribution  tt  the 
dynamical  matrix  where 

icfC)  =^^^ZOpfc;f  ( K) 

yields  the  full  dynamical  matrix  P^^-,KK^  : 

D^[]c\KK^)  =  D^[k\KK‘)+p^^  (16) 

The  normal  mode  frequ.encief  and  eigenvectors  were  obtained 
by  diagonalizing,  the  full  dynamical  matrix. 

i.  MODEL  PARAMETERS  AND  PHONON  DISPERSION  CURVES  FOR  PURE 
y-InSe 

-  ■  1 

The  central  potential  model  which  we  employ  contains 
twelve  parameters  to  be  determined:  A^,  .  — Ag,  .  i .  .Bg, 
|,*,Z2  .  The  experimental  data  that  are  available  for  y-InSe, 


riameiyv  infrared  and  Raman  spedtra,  are  insufficient  to- make 
a  ‘unidue  determination  of  these  parameters;.  We  iiave 
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Ffig.  2  Zone  center  normal  .mpdes  and  their  corresponding 
frequencies  for  -pure  y-InSe-i 


therefore  taken  as  a  starting  point  the  values  of  the 
correspohding  iparameters  for  GaS  obtained  by  Polian  [5]  by 
fitting  .phonon  dispersion  curves  (determineq  by  inelastic 
neutron  scattering)/,  elastic  cOnstant|>  and*  infrared  and 
Raman  data.  We  then  finestuneo  the  v|lues  of  the  parameters 
by  adjusting  them  to  reproduce  the  frequencies  and 
polarizations  of  the  optical  mopes  of  long  wavelength 
derived  from  infrared  and  ;Raman  data  fOr  yr^inSB  [jG]  which 
are  presented  in  FiCi  2.  The  values  of  the  model  parameters 
thus  determined  are  given  in  Ta_ble  1 . 


Table  1.  Short-range  and  Coulomb  parameters  for  y-InSe. 

The  units  for  the  A  and  B  parameters  are  N/m. 


i 

1 

2 

3 

4 

5 

Ai 

78.35 

32 . 17 

0.5 

7.0 

3.4 

Bi 

4.65 

31.0 

-0.196 

-8.584 

0.9 

★ 

0.71 

0.43 

— 

— 

— 

The  elements  of  the  dynamical  matrix  were  next 
evaluated  using  the  tabulated  values  of  the  model 
parameters.  Diagonalization  of  the  dynamical  matrix  for 
various  values  of  the  wave  vector  k  yielded  the  normal  mode 
frequencies  and  eigenvectors.  The  dispersion  curves  were 
constructed  for  several  high  symmetry  directions  in  the 
Brillouin  zone  and  are  presented  in  Fig.  3.  Since  there  are 
four  atoms  per  primitive  unit  cell  in  y-InSe,  there  ax  , 


A  r  B  r  z 

REDUCED  WAVE  VECTOR 


Fig.  3 


Phonon  dispersion  curves  for  pure  y-InSe. 


twelve  normal  modes  for  a  given  value  of  k  .  There  is 
considerable  dispersion  of  the  normal  modes  propagating  in 
directions  not  perpendicular  to  the  layers  as  a  consequence 
of  the  strong  intra-layer  binding  of  the  atoms.  However, 
the  weak  interlayer  binding  (note  the  small  values  of  and 
82  in  Table  1)  leads  to  nearly  flat  dispersion  curves  in  the 
direction  perpendicular  to  the  layers,  as  evident  in  Fig.  3 
for  the  rZ  direction. 

4.  MODEL  PARAMETERS  AND  PHONON  DISPERSION  CURVES  FOR 

LITHIUM  INTERCALATED  y-InSe 

The  intercalation  of  lithium  into  y-InSe  leads  to  new 

peaks  (6)  in  the  Raman  spectrum  at  92  cm~^  and  388  cm”^ .  In 

addition,  peaks  associated  with  pure  y-InSe  are  modified  by 

plasmon-phonon  coupling  (6) .  We  have  made  a  theoretical 

investigation  of  the  case  in  which  there  is  one  lithium  atom 

per  primitive  unit  cell  and  that  atom  occupies  a  site  of 

lowest  potential  energy  (A3-2  site  [7])  in  the  Van  der  Waals 

gap  between  layers .  We  assume  that  there  are  central 

potential  interactions  between  the  lithium  atom  and  the 

nearest  selenium  and  indium  atoms  with  interaction 

parameters  Ag,Bg  and  respectively.  Since  we  have 

only  two  new  normal  mode  frequencies  and  since  B  is 

typically  small  compared  to  A,  we  have  neglected  Bg  and  B.^ 

and  determined  Ag  and  h-j  by  associating  the  frequencies  92 
—  1  —1 

cm  and  388  cm  with  normal  modes  dominated  by  lithium 
atom  motion  parallel  and  perpendicular  to  the  layers, 
respectively.  The  resulting  values  for  Ag  and  A-y  are  1.61 
N/m  and  54.8  N/m,  respectively.  The  large  value  of  A^ 
relative  to  that  of  Ag  is  consistent  with  the  facts  that  the 
indium  atomic  radius  (0.76A)  is  larger  than  that  of  selenium 
(0.60A)  [8]  and  that  a  lithium  atom  in  the  A3-2  site  lies 
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directly  above  the  indium  atom  to  which  it  is  coupled.  The 
restoring  force  on  the  lithium  atom  is  therefore  dominated 
by  in  the  mode  in  which  the  lithium  atom  moves 
perpendicular  to  the  layers  {(0=388  cm~^)  and  by  Ag  in  the 
mode  in  which  the  lithium  atom  moves  parallel  to  the  layers 
(60=92  cm"^.) 

After  incorporating  the  force  constants  Ag  and  Aj  into 
the  dynamical  matrix,  we  calculated  the  normal  mode 
frequencies  and  eigenvectors  as  functions  of  the  wave  vector 
k  for  LiQ  g  InSe.  The  dispersion  curves  thus  obtained  are 
plotted  in  Fig.  4.  We  note  that  the  normal  modes  primarily 
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associated  with  the  InSe  lattice  are  little  affected  by  the 

presence  of  the  lithium,  since  the  mass  of  a  lithium  atom  is 

much  smaller  than  those  of  indium  and  selenium  atoms.  The 

normal  modes  primarily  associated  with  lithium  atom  motion 
—  1  —1 

at  92  cm  and  388  cm  are  essentially  flat  as  a  result  of 
the  weak  Li-Se  coupling  and  the  lack  of  coupling  between 
lithium  atoms  in  different  unit  cells . 

5,  CONCLUSIONS 

A  lattice  dynamical  model  has  been  developed  for  both 
pure  and  lithium-doped  y-InSe.  The  model  contains  short- 
range  and  Coulomb  interactions.  Phonon  dispersion  curves 
have  been  calculated  for  both  the  intercalated  and  non- 
intercalated  cases.  The  results  are  consistent  with 
observed  infrared  and  Raman  data. 
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D.  Molecular  Dynamics  Study  of  the  Lattice  Vibration 
Contribution  to  the  Frequency-Dependent  Dielectric  Constant 

of  Lithium  Iodide 


J.  Deppe*,  M.  Balkanski”,  R.  F.  Wallis  and  A.  R.  McGurn'®’ 
Physics  Department,  University  of  California,  Irvine 
Irvine,  California  92717 


Abstract 

A  molecular  dynamics  simulation  has  been  performed  on  the  crystal  lithium 
iodide,  Lil.  A  rigid  ion  potential  is  used  with  parameters  fit  to  thermal  ex¬ 
pansion,  isothermal  compressibility,  lattice  energy  and  the  frequency  of  the 
transverse  optical  mode  at  the  zone  center.  The  current-current  correlation 
function  has  been  calculated  at  T  =  200/^  and  400/^,  and  from  this  the 
absorption  and  dispersion  have  been  obtained.  .Anharmonic  broadening  is  ob¬ 
served  at  the  higher  temperature.  In  addition,  the  mean  square  displacements 
of  the  two  ions,  and  the  radial  distribution  function  are  calculated. 


1.  Introduction 


Lattice  dynamics  calculations  on  ionic  crystals  are  beset  by  a  number  of  interesting 
theoretical  problems.  The  most  evident  problems  w’hich  distinguish  the  study  of  ionic 
crystals  involve  the  proper  treatment  of  the  long  range  Coulomb  interactions  between 
ions.  This  difficulty  hcis  been  studied  over  a  period  of  many  decades  and  significant 
advances  in  the  handling  of  the  long  range  Couloumb  interaction  have  been  made  by 
EwaldL  Born*  and  other  workers  as  we  shall  detail  below.  In  addition  to  these  long 
range  electrostatic  effects,  atomic  polarization  effects  in  ionic  solids  have  been  treated 
by  Dick  and  Overhauser^,  Hardy"*,  and  Karo  and  Hardy®  using  shell  and  deformation 
dipole  models,  respectively.  The  successful  description  of  the  lattice  dynamics  of  the 
lithium  halides,  however,  is  still  hampered  by  problems  peculiar  to  these  particular 
ionic  compounds.  These  problems  include:  1)  the  large  non-central  interactions  which 
are  inferred  from  deviations  from  the  Cauchy  relation,  c\2  =  C44,  and  which  in  LiF^, 
has  the  value  of  cu  —  C44  =  —2.25  x  IQ'^^dynfcm?,  the  largest  of  all  the  alkali  halides; 
2)  the  importance  of  next  nearest  neighbor  short-range  interactions  due  to  the  small 
size  of  the  lithium  atom:  and  3)  the  large  anharmonic  contributions,  especially  for  the 
heavier  halide  salts.  In  an  attempt  to  address  a  solution  to  1)  through  3)  above,  V'erma 
and  Singh”  found  it  necessary  to  employ  a  modified  shell  model  which  included  3-body 
interactions®  to  calculate  the  phonon  dispersion  relations  of  the  lithium  halides.  .Along 
the  lines  of  Ref.  7,  Rastogi  et  al®  used  a  similar  model  to  calculate  the  phonon  dispersion 
relations  and  the  anharmonic  properties  of  the  lithium  halides.  Their  calculation  of  the 
real  and  imaginary  parts  of  the  self-energy  shift  at  u;  =  0  and  ij  =  ujto  indicates  that 
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the  anharniocic  contributions  increase  substantially  tvith  increasing  halide  size.  An 
important  conclusion  of  their  work  wais  that  the  perturbation  approach  employed  was 
legitimate  for  lithium  iodide  up  to  ISOiif,  which  we  believe  (see  below)  is  comparable 
to  or  less  than  Od^LU). 

In  this  paper,  we  investigate  the  temperature  dependence  of  the  infrared  lattice 
vibration  absorption  of  crystalline  lithium  iodide.  We  have  modified  a  rigid  ion  model 
first  proposed  by  Michielsen  et  al.'°  .Although  computer  resources  have  now  reached  the 
stage  such  that  more  sophisticated  potentials,  i.e.  3-body  interactions^  and  explicit 
polarization  terms*^,  can  be  employed,  the  use  of  a  central  force  interaction  potential 
is  justified  in  this  case  by  the  extremely  small  deviation  from  Cauchy  behavior  for 
lithium  iodide^^;  (ci2  —  +  c^^)  =  0.005.  In  addition,  although  polarization 

certainly  affects  the  dynamical  properties  of  this  crystal,  the  transverse  optical  mode 
at  fc  =  0,  which  is  the  mode  we  will  be  probing,  is  not  affected  to  a  large  extent. 
We  are  left  with  the  large  anharmonicity,  which  can  be  handled  nonperturbatively 
by  the  molecular  dynamics  method.  The  molecular  dynamics  simulation  is  suitable 
for  investigating  thermodynamic  and  dynamic  properties  of  materials.  It  has  been 
successfully  applied  to  fast  ion  conduction  and  defect  diffusion  in  crystalline  solids 
structural  investigations^®-*'  and  vibrational  studies*®  of  amorphous  solids,  and  studies 
of  correlated  motion  of  atoms  in  liquids*®.  One  of  the  main  reasons  the  use  of  this 
metho<i  has  become  so  widespread  is  the  ease  with  which  results  may  be  compared  with 
experimental  data.  In  principle,  any  experimentally  measurable  quantity  can  also  be 
extracted  from  the  molecular  dynamics  results.  In  the  present  treatment,  the  frequency 


and  temperature  dependent  conductivity  (which  is  measured  in  the  laboratory),  and 
from  this  the  dielectric  function,  is  directly  obtained  from  the  current-current  correlation 
function.  We  have  calculated  these  quantities  at  two  temperatures,  200K  which  we 
argue  is  on  the  order  of  the  Debye  temperature  for  this  material,  and  AQQK,  at  which  we 
see  large  enharmonic  broadening  in  the  correlation  functions.  Lowndes^®  has  measured 
the  transverse  optical  frequency  uJto  using  far  infrared  transmission  through  ultra  thin 
films  at  r  =  oK  and  T  =  Troom-  However,  we  have  not  uncovered  a  measurement  of 
the  temperature  dependent  broadening  with  which  we  can  compare  our  results.  Other 
properties  that  we  have  calculated  are  the  radial  distribution  function  and  mean  square 
displacements  of  the  ions. 

The  remainder  of  the  paper  follows  the  following  outline.  In  section  2,  the  crystal 
structure  and  relevant  properties  of  Lil  are  reviewed.  In  section  3  we  give  the  potential 
to  be  used  in  this  simulation,  a  variation  on  one  first  used  by  Michielsen  et  aP°.  The 
assumptions  and  constraints  inherent  in  the  model  will  be  enumerated.  The  equations 
of  motion  and  simulation  parameters,  such  as  the  timestep,  numerical  tolerances  and 
temperatures  investigated  are  also  given.  In  addition  the  determination  of  the  corre¬ 
lation  functions,  radial  distribution  function  and  mean  square  displacement  function  is 
discussed.  In  section  4,  results  are  presented  for  two  different  temperatures,  allowing 
an  estimation  of  the  magnitude  of  anharmonic  effects.  Finally,  in  the  last  section  are 
presented  conclusions  and  suggestions  for  future  study. 


2.  Properties  of  Lil 


Lithium  iodide  crystallizes  in  the  NaCl  structure.  The  bravais  lattice  is  face-centered 
cubic;  the  basis  consists  of  one  Li  atom  and  one  I  atom  separated  by  one-half  the  body 
diagonal  of  the  conventional  unit  cell.  Each  atom  haa  six  nearest  neighbors  of  the  op¬ 
posite  species.  The  reciprocal  lattice  vectors  for  an  fee  lattice  are  given  by  the  primitive 
translation  vectors  of  a  body-centered  cubic  lattice.  The  nearest  neighbor  distance  is 
r„„  =  3.OOA  at  room  temperature^*,  and  the  cube  edge  has  length  a  =  2r„„.  There  is  a 
relative  lack  of  experimental  and  theoretical  work  on  this  crystal;  the  former  is  possibly 
due  to  the  fact  that  the  solid  is  extremely  hygroscopic.  We  have  uncovered  experimental 
data  on  the  melting  temperature*®,  thermal  expansion  and  isothermal  compressibility^^, 
the  high  frequency  dielectric  constant^®,  and  lattice  energy*®.  The  transverse  optical 
frequency  at  the  zone  center^  has  been  determined  using  far  infrared  transmission  mea¬ 
surements  through  ultra  thin  films;  however,  the  variation  of  the  absorption  peak  with 
temperature  heis  not  been  discussed  to  our  knowledge.  A  measurement  of  the  elastic 
constants*^  has  yielded  cu  «  c^^,  indicating  that  a  central  force  approximation  may  be 
a  very  good  one.  The  other  lithium  halides,  in  contrast,  show  large  deviations  between 
C12  and  the  difference  for  LiF  is  the  largest  of  all  the  alkali  halides.  There  has  also 
been  a  measurement  of  the  ionic  conductivity  of  LilP  The  study  of  ionic  conductivity 
in  the  alkali  halides  has  a  long  history  and  is  known  to  occur  almost  exclusively  through 
the  diffusion  of  positive  ion  vacancies^^.  These  vacancies  arise  either  through  thermal 
activation,  leading  to  pairs  of  positive  and  negative  ion  vacancies  known  ais  Schottky 
pairs,  or  as  a  result  of  doping.  In  our  case,  there  are  no  defects  due  to  doping,  and 
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TABLE  I:  Value  of  the  Debye  temperature  for  some  alkali  iodide  compounds. 


Substance 

0D 

.Nal  1 

164K 

KI 

132K 

Rbl 

103K 

diffusion  would  be  due  entirely  to  Schottky  pairs.  However,  it  can  be  expected  that, 
although  ionic  diffusion  will  be  quite  small  in  the  solid  phase  of  Lil  without  the  intro¬ 
duction  of  impurities  or  defects,  at  temperatures  not  too  far  below  the  melting  point 
the  lithium  may  indeed  be  able  to  diffuse  through  interstitial  hopping.  The  iodine  ion  is 
huge  in  comparison  to  the  lithium  ion,  and  may,  in  effect,  leave  large  channels  through 
which  the  lithium  can  move.  In  the  former  case  of  diffusion  by  Schottky  defects,  the 
mobility  activation  energy  and  defect  formation  activation  energy  have  been  measured 
by  Jackson  and  Young*^  to  be  0.43  and  l.OOeV  respectively. 

The  crystal  is  single  phased  throughout  the  temperature  range  up  to  its  melting 
point,  which  is  1A2K.  .Although  we  have  been  unable  to  find  a  value  for  its  Debye 
temperature,  it  can  perhaps  be  estimated  by  considering  the  value  of  this  quantity  for 
the  other  alkali  halides.  In  Table  I,  we  list  the  Debye  temperature  for  three  alkali- 
iodides*®’^®.  By  simply  extrapolating  this  data  to  lithium  iodide,  we  might  therefore 
expect  that  the  value  of  Od  is  of  the  order  of  200K.  .An  analysis  which  takes  into 
account  the  variation  in  sound  speed  and  lattice  spacing  indicates  that  6D{LiI)  may  be 
somewhat  lower,  on  the  order  of  ISOA'. 


3.  Theory 

3.1  Interatomic  potential 


In  this  work,  we  have  chosen  a  modified  version  of  an  interatomic  potential  proposed 
by  Michielsen  et  al*°.  The  largest  term  in  the  potential,  in  terms  of  its  contribution  to 
the  ground  state  energy,  is  unambiguously  defined,  since  Lil  is  an  ionic  crystal.  The 
Coulomb  energy  for  this  solid  is  found  to  account  for  «  90%  of  the  total  potential  energy, 
as  is  typical  for  all  the  alkali  halides.  A  modified  Bom-Mayer  term  is  used  for  the  short 
range  repulsive  term.  In  addition,  dipole-dipole  and  quadrupole-dipole  attractive  terms 
and  a  hard  core  repulsive  term  are  included.  The  pair  potential  is  given  by, 

•ine(^r{j  r--  Vij  rfj  rf- 


(1) 


where 

r,-j  =  separation  of  atoms  i  and  j, 
Zi  =  valence  of  atom  i. 


Tyo  =sum  of  ionic  radii^", 

dj  =coefficient  of  dipole-dipole  interaction, 

dij  =coefficient  of  dipoie-quadnipole  interaction, 

=coetficient  of  hard  core  repulsion, 

6,  k  =parameters  fit  to  thermal  expansion  and  isothermal  compressibility. 
We  can  express  the  total  potential  energy  V’(r)  as 
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where  A  is  the  Madelung  constant,  N  is  the  number  of  -*.it  cells,  r  is  the  nearest 
neighbor  distance,  and  it  is  understood  that  the  sum  does  not  include  the  term  i  =  j. 
We  have  used  the  method  of  Sangster  and  Dixon^®  in  the  evaluation  of  the  long-range 
Coulomb  term  with  the  Ewald  summation  technique. 

There  are  two  parameters  in  this  potential,  both  in  the  modified  Born-Mayer  term, 
which  are  fit  to  experimental  data:  the  prefactor  6,  and  the  inverse  decay  length  k.  We 
may  determine  these  parameters  by  considering  the  following  two  expressions^®  which 
relate  derivatives  of  the  potential  at  non-zero  temperature  to  thermodynamic  quantities. 


,f±. 

dr'^ 


2“  I  _  9Vm  f-i  f  f  \  I  Toifd^^  1.  ^T/v^ 


(3) 

(4) 


where  $  is  the  energy  density,  tq  is  the  equilibrium  lattice  constant,  a  is  the  thermal 
expansion  coefficient,  is  the  isothermal  compressibility  and  V\f  is  the  volume  per 
molecule. 

The  following  are  experimental  values  for  the  lattice  energy  U,  thermal  expansion 
coefficient,  the  isothermal  compressibility  and  two  of  it’s  derivatives. 


-U  =  754.3)fcJ,‘° 
a  =  16.7  X 
jS  =  5.83  X 

^{^)t  =  37.3  X  lO-V/N,^ 


(5) 


8 
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wpo  =  2.7  X 

We  desire  that  the  potential  energy  satisfy  the  following  objectives:  reproduce  the 
lattice  energy  of  the  crystal;  reproduce  the  frequency  of  the  TO  mode  at  the  zone  center; 
lead  to  the  satisfaction  of  Eqs.  (3)  and  (4).  In  addition,  we  of  course  require  that  the 
ions  be  repulsed  as  they  approach  each  other,  and  this  sets  rather  stringent  limits  on 
the  hard  core  term  Cijlrfj. 

We  fit  the  quantities  listed  above  with  the  following  parameters:  the  magnitude  of  . 
the  prefactor  to  the  Pauli  repulsion,  b;  the  decay  constant  k;  and  the  magnitudes  of 
the  dipole  interactions.  The  strength  of  the  hard  core  repulsion  between  the  atoms 
e,j  is  in  effect  fixed.  This  Izist  statement  is  true  because  the  repulsion  must  be  strong 
enough  to  prevent  a  double  well  from  developing  between  ions  (these  arising  from  the 
Coulomb  attraction  and  possibly  the  dipole  interactions),  and  must  be  small  enough  to 
keep  the  hard  core  deep.  We  have  chosen  to  vary  the  strength  of  the  dipole  interactions 
from  the  values  given  by  Ladd^'  because  the  frequency  of  the  TO  mode  at  9  =  0 
is  considerably  less  than  what  this  model  predicts,  and  the  very  large  coefficients  of 
the  dipole  interactions  are  the  only  parameters  which  can  be  varied  without  adversely 
affecting  the  lattice  energy  and/or  the  two  thermodynamic  relations,  Eqs.  (3)  and  (4). 

The  procedure  is  the  following;  first  the  lattice  spacing  is  set  to  3i4  at  room  tem¬ 
perature,  then  we  adjusted  the  b  and  k  to  satisfy  the  thermodynamic  relations  given 
in  Eqs.  (3)  and  (4).  It  should  be  noted  that  these  relations  are  not  very  sensitive  to 
the  c,j  and  dij.  The  TO  frequency  is  then  calculated,  we  then  adjust  the  Cij  and 

-j 
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TABLE  II:  Values  of  the  parameters  used  in  Lil  calculations,  m  addition,  i:  =  1.1  x  10^°m“^,6  = 
2.8  X  10-58jm<  _ _ 


Li-l 

Li-Li 

i-I 

nio(A) 

3.00 

0.95 

2.05" 

“c.v7l0-^(j-m^l 

25.0 

8.54 

266.0 

-  m«) 

114. 

36.8 

1000.0 

e,j7lO~^‘”(J  -  m^) 

6.5 

10.0 

400.0 

re-determine  b  and  k,  following  the  same  steps  until  the  TO  frequency  is  close  to  the 
experimental  value  2.7  x  10‘^rad— s"'.  Finally,  the  e,y  are  adjusted  to  ensure  that  there 
is  no  attraction  between  the  ions  as  they  approach  closer  than  3A. 

We  find  the  lattice  energy  to  he  U  —  725kJ,  within  5%  of  the  experimental  value 
of  IbikJ,  Eq.  (5).  In  addition  the  calculated  TO  frequency  is  %  3.1  x  lO'^rad  — 
within  10%  of  the  experimental  value  of  2.7  x  lO^^rad  —  5“^  The  expressions  involving 
the  first  and  second  derivatives  of  the  potential  are  fit  exactly  at  room  temperature, 
and  within  5%  error  at  all  others.  We  have  varied  only  the  lattice  constant  as  a  function 
of  temperature,  leaving  one  parameter  and  two  equations  to  be  fit.  This  has  been  done 
by  splitting  the  difference  between  the  results  from  the  two  equations. 

We  have  found  that  restricting  the  short  range  interactions  to  the  second  nearest 
neighbors,  as  has  been  done  often  in  the  past,  is  insufficient  here.  Neglecting  terms 
beyond  second  neighbors  was  found  to  result  in  a  5%  discrepancy  compared  to  the 
result  including  their  contribution.  Therefore,  all  atoms  out  to  the  tenth  neighbors  are 
included  in  this  part  of  the  calculation.  The  values  of  the  parameters  used  are  listed  in 
Table  II. 


There  are  certain  defects  inherent  in  such  a  model.  Among  them  are: 


1)  Tunneling  effects,  which  may  well  be  important  considering  how  small  the  lithium 
ion  is,  and  the  importance  of  interstitial  diffusion  in  this  system,  are  neglected. 

2)  The  potential  is  that  of  a  rigid  ion  model,  so  there  are  no  polarization  effects 
present.  We  are  interested  in  the  response  at  the  transverse  optical  (TO)  frequency, 
however,  which  is  not  affected  to  the  same  extent  that  the  longitudinal  optical  (LO) 
mode  is.  .Although  simulations  have  been  successfully  performed  using  rigid  ion  models 
for  extremely  polarizable  materials  in  the  past^^’“,  {Agl),  it  would  be  very  desirable  to 
do  a  molecular  dynamics  simulation  using  a  shell  model  potential  on  this  system. 

3)  Since  the  parameters  of  the  potential  are  fit  to  low  temperature  experimental 
data  for  which  the  harmonic  approximation  gives  a  complete  description  of  the  system 
dynamics,  the  potential  is  only  accurate  near  the  minimum  of  the  potential  energy;  we 
therefore  might  expect  problems  at  very  high  temperatures,  i.e.  in  the  liquid  phase. 

We  note  here  that  the  original  potential  used  in  Ref.  10  did  not  have  the  hard  wall 
repulsion  and  therefore  did  not  prohibit  pairs  of  ions  from  collapsing  on  themselves.  This 
difficulty  is  only  critical  at  very  high  temperatures,  where  the  model  crystal  is  entirely 
unstable.  The  absence  of  a  repulsive  term  in  the  potential  that  can  counterbalance  the 
Van  der  Waals  terms  at  small  separations  leads  to  a  collapse  of  the  pair. 

3.2  Equations  of  Motion 

We  treat  a  crystal  supercell  in  the  form  of  a  cube  containing  216  atoms.  Periodic 
boundary  conditions  are  imposed.  The  classical  system  of  equations  to  be  solved  <'on- 
sists  of  216  X  3  =  648  coupled  differential  equations,  3  (corresponding  to  the  cartesian 
directions)  for  each  of  the  216  atoms.  The  long-range  Coulomb  force  couples  atoms  at 
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all  distances  and  is  handled  by  the  Ewald^  method.  The  reciprocal  space  summation 
effectively  accounts  for  the  interactions  of  atoms  which  are  separated  by  more  than  one 
supercell  distance,  which  in  our  case  is  IsA.  We  have  solved  this  system  of  equations, 
which  include  anharmonic  contributions  exactly,  to  a  numerical  accuracy  of  1  part  in 
10^ 

We  begin  with  Newton’s  equation,  a  second  order  differential  equation. 

Fa  =  mxa,  (6) 

and  express  this  as  two  first  order  coupled  ordinary  differential  equations,  (we  now  drop 
the  cartesian  subscript  for  convenience) 

p=%  (V 

P  =  n,§.  (8) 

The  force  arises  from  interactions  involving  all  atoms  in  the  supercell  and  beyond,  and 
is  in  addition  a  vector  quantity.  We  indicate  a  vector  quantity  with  boldface. 

Now  let  us  make  the  following  change  in  scale: 

m  =  m'mo,  (9) 

X  =  x'xo, 

F  =  F'Fo, 

where  mo  =  mp  is  the  proton  mass,  lo  =  and  Fo  =  We  also  write 

Cij  =  cO-  X  10-"9(J  -  m®),  dij  =  4  X  -  m®),  e.y  =  e?-  x  -  m®)  and 

b  =  b°  X  —  m*).  Using  p'  =  p/>/moXoFo.  and  t'  =  Faj xoniQ,  the  equations 
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of  motion  can  be  expressed  in  terms  of  quantities  which  are  all  of  order  unity,  which  is 
necessary  from  a  computational  standpoint.  We  have 


p'  -  ‘‘p' 

^  -TF' 


(10) 

(11) 


If  these  scales  are  combined,  we  find  that  the  "real”  time  is  related  to  the  scaled  time  by 
t'  =  2.42  X  We  have  employed  a  scaled  timestep  of  1  in  all  simulations  presented 
here,  which  corresponds  to  4.1  femtoseconds. 

Now  that  the  time  scales  are  established,  it  will  be  necessary  to  scale  the  velocities, 
or  momenta,  in  order  to  set  the  effective  temperature  of  the  crystal.  We  have  that 
Ek  =  \NkBTtff,  where 


Ek  =Ep?/(2"‘.)- 

t 

In  terms  of  the  scaled,  or  primed  variables,  this  becomes. 


Ek  = 


^(moFoXo)pO^ 

2m\mo 


=  io-‘8y"i£0! 
r  2m<  ’ 


We  can  solve  for  Te//,  using  ks  =  1.38  x  and  find 

2Ek 


ZNke 


=  223.6£;;f. 


(12) 


(13) 


(14) 
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3.3  Connections  with  experimental  data 

VVe  now  turn  attention  to  the  current-current  correlation  function,  which  will  be  shown 
to  be  related  to  the  conductivity  and  dielectric  constant.  First,  we  define  the  frequency 
dependent  current  density-current  density  correlation  function  as 

na(w)  =  f  dte''^^{ja{t)ja{0)),  (15) 

J  —<X> 

where  ia{t)  is  the  g  =  0  fourier  component  of  the  current  density  vector  and  the  angular 
brackets  denote  an  ensemble  average.  We  will  refer  to  the  current  density-current 
density  correlation  function  as  simply  the  current-current  correlation  function.  We  can 
express  this  as^"* 

n.(a-)  =  Um  i  J' 

These  integrations  cannot  be  performed  exactly  as  they  extend  to  co;  however,  cutting 
off  the  summation  at  t  =  r  will  give  the  resulting  correlation  function  a  width,  Aw 
=»  27r/r,  which  can  be  compared  with  the  actual  width. 

Many  advanced  texts  on  many-bodj'  physics  discuss  the  connection  between  this 
correlation  function  and  the  frequency  dependent  conductivity.  The  real  part  of  the 
conductivity  at  9  =  0  can  be  expressed  in  terms  of  the  velocity-velocity  correlation 
function  as^® 

Re[aay{ij)]  =  2^^(1  -  J  ^dte''^\va{t)v^{Q)) a,  (16) 

where  the  subscript  G  indicates  an  averageover  a  grand  canonical  ensemble,  0  =  IfkgT, 
Q  is  the  particle  charge,  V  is  the  crystal  volume  and  h  =Planck’s  constant  divided  by 
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2n.  Presently,  we  are  onlj’  considering  the  diagonal  component  of  this  tensor,  we  have 
therefore  used  the  averar/- t*  three  cartesian  components  to  improve  our  statistics. 
All  results  presented  ' .  ’j*  -u  averaged  in  this  way.  We  can  move  the  factors  of 

Q  and  V  into  the  grand  can  .acal  ensemble  average  to  obtain  the  conductivity  in  terms 
of  the  current  density-curr ,  v  'density  correlation  function.  Comparing  Eqs.  (15)  and 
(16),  we  find  that  the  fur.- : ,  :  '.ve  are  calculating  is  related  to  the  frequency  dependent 
conductivity  by 

n(cc;),(l  -  e-^«“)F/(2^u;)  =  Re[a^,,{uj)].  (17) 

We  can  also  relate  these  expressions  to  the  real  and  imaginary  parts  of  the  dielectric 
function,  which  are  more  amenable  to  experimental  determination. 

We  have 

e(a.)  =  (£„  +  i— );  (18) 

UJ 

therefore, 

/m(e(a;))  :=  2Tn(a;)(l  -  (19) 

We  can  then  find  the  real  part  of  £(a’)  using  the  Kramers-Kronig  relations. 

One  final  observation  concerning  the  calculation  of  the  current-current  correlation 
function  is  the  following:  Many  quantities  derived  from  molecular  dynamics  simulations 
have  within  them  some  fluctuations,  either  from  small  temperature  variation: ,  round-off 
error,  or  a  lack  of  suitably  long  time  averaging.  This  effect  is  particularly  pronounced  in 
correlation  functions,  which  sample  over  long  distances  in  the  crystal.  It  has  therefore 
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become  commonplace^,  if  not  standard,  to  smooth  correlation  functions  with  a  gaussian 
w  hich  is  centered  on  each  data  point  and  whose  width  is  nominally  some  small  quantity. 
In  the  results  to  follow,  we  have  employed  such  a  smoothing  procedure  using  the  function 

r  (20) 

y/TTO  oo 

where  wo  is  the  frequency  of  interest  and  6  =  o."'  <  lO^^rad  —  s“^  which  is  3/500  of  the 
full  frequency  width  over  which  the  correlation  i  anctions  are  calculated. 

Next  we  define  the  mean  square  displacement  as  a  function  of  time,  MSD{t),  which 
is  related  to  the  self  diffusion  of  the  two  ion  types  in  this  case, 

MSD{t)=^{{ri{to  +  t)-ri{to))^).  (21) 

The  angular  brackets  indicate  an  average  over  all  atoms  or  ions  of  type  i. 

The  direct  connection  between  the  mean  square  displacement  and  the  diffusion  co- 
efBcient  D  for  that  particular  particle  is  given  by^^, 

=  6Di,  (22) 

where  t  is  the  time. 

We  also  calculate  the  radial  distribution  function,  RDF.  This  quantity  simply  ex¬ 
presses  the  average  number  of  atoms  within  some  distance,  r,  of  some  particular  atom. 

4.  Molecular  Dynamics  Simulation 

In  this  section  we  give  results  from  molecular  dynamics  calculations  on  Lil  at  two 
temperatures.  We  have  calculated  the  current-current  correlation  function,  the  mean 
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square  displacements  for  the  two  atoms,  and  the  radial  distribution  function.  The 
current-current  correlation  function  is  related  to  the  ionic  conductivity  through  a  factor 
of  (1  —  We  then  use  the  relation  between  the  ionic  conductivity  and  the 

imaginary  part  of  the  dielectric  function  to  find  the  absorption.  The  real,  or  dispersive 
part  of  the  dielectric  function  is  then  calculated  from  the  imaginary  part  using  the 
Kramers-Kronig  relations.  We  have  obtained  these  quantities  at  the  temperatures  200/if 
and  400/f ,  the  first  being  well  below  the  melting  temperature  Tm,  which  is  742/f ,  and 
of  the  same  order  cf  magnitude  as  the  Debye  temperature(0£))"  The  second  temperature 
is  still  well  below  hut  presumably  almost  29d  and  should  exhibit  large  anharmonic 
effects.  We  account  for  the  thermal  expansion  of  the  crystal  using  Eqs.  (3)  and  (4). 

The  thermal  expansion  produces  a  decrease  in  the  TO  frequency  cis  a  function  of 
increasing  temperature.  This  is  simple  understood  as  a  consequence  of  the  atoms 
moving  further  from  the  absolute  minimum  of  their  potential  energy  positions  bls  the 
temperature  increases.  In  addition  to  this  "static”  frequency  shift,  there  is  also  to  be 
expected  a  dynamic  downward  shift  in  frequency  with  increasing  temperature  due  to 
the  vibrating  atoms  reaching  into  regions  of  the  potential  with  decreasing  curvature. 

It  is  not  at  all  clear  that  the  initial  positions  and  momenta  represent  an  equilibrium 
configuration;  therefore,  we  integrate  the  equations  of  motion  ahead  for  some  thermal- 
ization  time,  Ith-  One  criterion  for  determining  when  this  thermalization  time  has  been 
reached,  and  the  one  to  be  used  here,  is  that  the  effective  temperature  Te/j  =  EkH^ks), 
varies  by  less  than  a  couple  of  percent  across  the  sample,  and  that  it  not  vary  with  time. 
We  have  tricked  the  system  into  quickly  achieving  this  thermal  equilibrium  in  the  fol- 
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lowing  way;  at  the  beginning  of  the  simulation,  the  potential  and  kinetic  energies  are 
calculated,  the  momenta  are  scaled  to  correspond  to  the  desired  temperature,  tnen  the 
positions  are  adjusted  so  that  the  kinetic  and  potential  energies  have  roughly  the  same 
value,  as  they  should  according  to  the  equipartition  theorem.  Although  this  does  not 
result  in  perfect  thermal  equilibrium,  it  decreases  the  computer  time  necessary  to  reach 
this  stage.  We  used  18,000  timesteps  to  reach  a  good  thermal  equilibrium  at  200if  and 
400/C. 

4.1  Results  and  Discussion 

We  performed  simulations  on  a  216  atom  cube  with  periodic  boundary  conditions.  The 
cutoff  parameter  in  the  Madelung  sum  wcis  set  to  0.35,  and  we  sum  over  reciprocal  lattice 
vectors  up  to  relative  magnitude  =  48.  The  energy  and  temperature  fluctuated  less 
than  4%  over  the  entire  simulation  time  for  all  the  Ccises  presented  here.  The  timestep 
was  4.1  femtoseconds.  This  timestep  is  somewhat  arbitrary  as  the  differential  equation 
solver  uses  as  many  timesteps  as  is  necessary  to  meet  the  tolerance  requirements  as  set 
by  the  user,  which  Wcis  set  to  1  part  in  10^.  We  have  performed  an  ensemble  average 
over  two  sets  of  different,  but  random,  initial  conditions.  All  of  the  results  presented 
below  have  been  averaged  over  these  two  separate  simulations. 

We  first  present  results  for  T  =  200/C,  for  which  18000  timesteps  were  required 
for  thermalization.  The  lattice  spacing  and  TO  frequency  at  the  zone  center  for  this 
temperature  are  2.99A  and  3.4  x  lO'^rad  —  s“*,  respectively.  The  correlation  func¬ 
tions  were  calculated  over  8000  timesteps  and  averaged  over  two  sets  of  random  initial 
conditions.  These  numbers  correspond  to  73.8  x  10"*^  seconds  of  thermalization  time 


and  32.8  x  10“^^  seconds  of  data  gathering  time.  This  is  indeed  a  rather  short  time, 
however,  the  correlation  functions  show  little  or  no  change  after  6000  steps.  We  first 
present  in  Fig.  1  the  imaginary  part  of  the  dielectric  constant,  or  the  absorption  peak 
from  the  current-current  correlation  function  at  T  =  200if ,  in  arbitrary  units.  Note 
the  sniall  oscillations  superimposed  on  the  main  peak;  these  could  be  smoothed  by 
widening  the  width  of  the  smoothing  gaussian  as  described  in  the  previous  section.  We 
will  encounter  these  features  in  the  data  for  T  =  400K  as  well,  and  will  take  up  their 
possible  causes  and  solutions  then.  The  function  peaks  at  the  TO  frequency  at  9  =  0, 
3.4  X  10‘^rad  —  as  expected,  and  shows  a  small  tail,  ais  seen  in  Fig.  1. 

Using  the  Kramers-Kronig  relations^,  we  have  calculated  the  dispersion,  or  real  part 
of  the  dielectric  function,  with  the  results  for  8000  timesteps  at  T  =  200/^  also  shown 
in  Fig.  1.  As  is  typical,  the  region  of  anomalous  dispersion  coincides  with  the  region  of 
strong  absorption. 

In  Fig.  2  is  plotted  the  radial  distribution  function  for  1500  timesteps.  The  results 
for  longer  run  times  and  for  T  =  400iif  are  indistinguishable.  It  is  interesting  to 
note  the  extremely  sharp  and  well  defined  peaks  correspKjnding  to  the  various  neighbor 
positions.  In  Table  III  we  give  the  number  of  neighbors  and  their  distances  for  the 
rock-salt  structure  with  the  nearest  neighbor  distance  normalized  to  unity.  We  can 
identify  up  to  the  10‘*  neighbor  from  the  graph  and  see  that  the  area  under  each  peak 
reproduces  the  number  of  the  various  neighbors  very  precisely. 

Next,  in  Fig.  3  are  given  the  mean  square  displacements  for  Li'*'  and  I~  ions  at  200K. 
We  note  immediately  that  there  is  no  diffusion,  ris  is  expected  at  this  low  temperature. 
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Figure  i*  Real  and  imaginary  parts  of  the  dielectric  function  of  Lil  (in  arbitrary  units)  calculated  over 
8000  timesteps.  This  data  represents  an  average  over  two  sets  of  random  initial  conditions  calculated 
at  a  temperature  of  200K.  The  lattice  spacing  and  TO  frequency  at  the  rone  center  are  2.99A  and 
3.4  X  iO^^rod  — respectively. 
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Figure  3:  Mean  square  displacement  as  a  function  of  time  at  T  =  200A'  for  iodine,  solid  line,  and 
lithium  dashed  line. 
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TABLE  III:  Numbers  of  various  neighbors  and  their  distances  for  Lil 


Neighbor  number 

number  of  neighbors 

normalized  distance 

1 

6 

1 

2 

12 

y/2  «  1.41 

3 

8 

v/3  «  1.73 

4 

6 

2 

5 

24 

6 

24 

7 

12 

v/8  «  2.83 

8 

27 

3 

9 

12 

■BSfEBSSlil 

10 

12 

If  self  diffusion  were  taking  place,  we  would  see  the  mean  square  displacement  increase 
linearly  with  time,  as  indicated  in  Eq.  (22).  We  can  estimate  the  frequency  of  oscillation 
of  the  lithium  using  the  spacing  of  the  peaks,  which  is  150  femtoseconds  for  the  lithium 
and  600  femtoseconds  for  the  iodine.  The  frequency  of  the  iodine  motion  is  close  to 
one-quarter  that  of  the  lithium,  which  corresponds  very  closely  to  the  square  root  of 
the  ratio  of  their  masses,  which  is  \Jmulmi  =  0.23. 

We  now  turn  attention  to  a  simulation  on  the  same  system  but  at  AOOK.  One 
advantage  of  performing  these  calculations  at  a  temperature  larger  than  the  Debye 
temperature  is  that  any  inhomogeneities  in  the  initial  conditions  should  be  evened 
out  quickly  due  to  the  large  anharmonicities  at  work.  We  have  used  the  same  initial 
positions  and  momenta  as  in  the  previous  simulation,  but  scaled  to  correspond  to  this 
temperature.  The  interatomic  spacing  and  TO  frequency  at  400/<'  are  ss  3.02A  and 
2.99  X  lO'^rad  —  respectively.  18,000  timesteps  were  used  for  thermalization,  and 
the  same  timestep  of  4.1  femtoseconds  v/as  used  throughout. 
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In  Fig.  4  are  plotted  the  absorption  and  dispersion  after  8000  timesteps  at  400i£r. 
The  gross  structure  of  these  plots  exhibits  the  following  features.  The  middle  of  the 
three  maxima  occurs  very  near  the  TO  frequency  of  2.99  x  lO^^rad  —  although 
its  exact  position  cannot  be  pinned  down.  The  width  of  the  function  has  increaised 
noticeably,  from  a  full  width  at  half  maximum  of  «  0.2  x  lO'^rad  —  at  T  =  200ii£!’, 
in  Fig.  1,  to  «  1.2  X  lO^^rad  —  s~^  in  Fig.  4.  In  addition,  the  length  of  the  high 
frequency  tail  is  greatly  increased  at  the  higher  temperature. 

One  disturbing  feature  of  these  correlation  functions  is  the  persistence  of  the  periodic 
maxima  and  minima.  It  is  found  that  temporal  and  ensemble  averaging  decreases  the 
magnitude  of  these  oscillations.  Therefore,  we  believe  this  small  scale  structure  to  be 
unphysical  and  a  consequence  of  the  finite  time  averaging. 

Finally  in  Fig.  5  is  plotted  the  mean  square  displacement  calculated  after  8000 
timesteps.  The  magnitude  of  the  displacements  for  the  two  species  is  the  same  indicating 
that  the  temperatures  of  the  two  sublattices  are  the  same. 

5.  Conclusion 

We  have  calculated  the  frequency  dependent  dielectric  function,  the  radial  distribution 
function  and  mean  square  displacements,  all  as  functions  of  temperature,  for  crystalline 
lithium  iodide.  A  rigid  ion  model  Weis  used  for  the  interatomic  potential.  The  various 
parameters  in  this  potential  were  fit  to  the  thermal  expansion,  the  isothermal  compress¬ 
ibility,  the  lattice  energy  and  the  frequency  of  the  TO  mode  at  the  zone  center.  We 
have  also  included  a  term  which  acts  as  a  hard  core  repulsion,  preventing  the  solid  from 
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Figure  4:  Real  and  imaginary  parts  of  the  dielectric  function  at  T  =  AQOK  after  8000  timesteps  of 
data  taking.  The  lattice  spacing  at  this  ;c.-ij>erature  is  3.02A.  This  absorption  peaks  near  the  TO 
frequency  of  2.99  x  lO^^rad-  s"^  however,  the  broadness  inhibits  an  exact  determination  of  the  peak 
position. 
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Figure  5:  Mean  square  displacement  as  a  function  of  time  at  T  =  400/f  for  iodine,  solid  line,  and 
lithium  dashed  line. 
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collapsing  at  high  temperatures. 

A  Icirge  increase  in  anharmonic  broadening  in  the  frequency  dependent  correlation 
functions  was  observed  in  passing  from  T  =  200/^  to  T  =  400/^,  indicating  that  the 
Debye  temperature  for  this  material  most  likely  lies  in  this  range,  as  we  had  inferred 
from  data  from  other  alkali  iodides  (see  Table  I).  We  have  been  able  to  reproduce 
the  TO  frequency  within  10%  of  the  experimental  value,  although  no  experimental 
measurement  of  the  temperature  dependence  of  the  infrared  absorption  or  dispersion 
for  this  material  has  been  found.  It  would  be  very  interesting  if  experimental  data  on 
these  temperature  dependences  were  to  become  available  to  provide  a  comparison  with 
our  theoreticcil  results. 
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